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Ab  initio  self-consistent  field  (SCF)  calculations  of  the  atonic 
polar  tensors  (APTs)  for  all  the  atoms  of  the  HOX  series  of  triatomic 
molecules  have  been  carried  cut  with  the  GAUSSIAN  70  program  using  a 
4-31G  basis  set  of  wavefunctions.  Here  the  X has  been  Ci , 0,  H or  F. 
We  have  compared  these  ab  initio  APTs  with  those  obtained  by  trans- 
ferring the  APTs  for  these  atoms  from  experimental  values  measured 
for  other  molecules. 

The  ab  initio  and  transferred  APTs  have  been  used  to  predict  the 
infrared  intensities  for  all  of  the  fundamental  bands  of  the  HOX 
series.  The  theoretical  spectra,  predicted  using  ab  initio  APTs 
and  also  frcm  transferred  APTs,  have  been  compared  with  the  experi- 
mental spectra  for  each  of  these  molecules. 

Our  intensity  calculations  based  on  the  ab  initio  APTs  are  on 
the  same  level  of  accuracy  as  those  obtained  using  transferred  APTs. 
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We  have  found  that  the  intensities  obtained  either  from  ab  initio 


SCF  calculations  of  APTs  or  from  transferred  APTs  are  in  good  agree- 
ment with  the  experimental  values  for  the  absolute  intensities, 
where  available,  or  with  the  relative  intensities  provided  the 
APTs  are  properly  transferred  from  atoms  in  similar  chemical  environ- 
ments. 

In  addition,  we  have  calculated  the  charge,  charge-flux  and 
overlap  contributions  to  the  ab  initi o APTs  for  all  of  the  atoms  of 
the  HOX  molecules.  A correlation  has  been  found  between  the 
corresponding  elements  of  the  charge-flux  tensors  for  the  hydrogen 
atom  of  the  HOX  series  and  the  Mulliken  electronegativity  of  the  X 
atom. 

Absolute  intensities  calculated  for  HOC#,  both  from  the  ab  initio 
SCF  calculation  of  APTs  and  from  the  transferred  APTs  are  compared 
in  some  detail  with  available  experimental  data.  For  HO^  and  HOF 
the  relative  intensities  frcm  the  calculations  have  been  compared 
with  the  relative  experimental  intensities  frcm  matrix  spectroscopic 
studies.  The  latter  are  fit  by  a computer-simulated  spectrum  obtained 
as  a combination  of  Gaussian  functions. 

Our  intensity  calculations  for  HOCS,  and  HO^  are  expected  to  be 
very  useful  for  the  application  to  spectroscopic  infrared  measurement 
of  the  concentration  of  those  species  in  the  earth's  stratosphere. 
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CHAPTER  1 
INTRODUCTION 

Since  the  theory  of  atomic  polar  tensors  (APTs)  was  introduced 
by  Biarge,  Herranz  and  Morcillo  [1]  and  rediscovered  and  extended 
by  Person  and  Newton  [2] , the  chapter  on  studies  of  intensities 
in  infrared  spectroscopy  has  taken  a new  dimension  with  considerable 
promise  for  the  future  in  this  field.  It  is  sufficient  to  say 
that  those  ingenious  and  revolutionary  ideas  make  it  now  possible 
to  predict  successfully  the  infrared  intensities  of  a number  of 
molecular  species,  using  intensity  parameters  which  can  be  trans- 
ferred from  chemically  related  molecules,  aided  by  quantum  mechan- 
ical calculations. 

This  idea  has  always  been  much  alive  in  the  minds  of  those 
who  have  dreamed  to  analyze  the  infrared  spectrum  of  a molecule 
utilizing  a powerful  framework  of  theory  which  enables  one  not  only 
to  predict  thek  experimental  frequencies  but  also  to  calculate  the 
infrared  intensities  for  each  fundamental  mode  of  the  molecule. 

The  beauty  of  the  theory  of  APTs  is  in  the  fact  that  only  very 
few  parameters  are  used  in  the  prediction  of  infrared  intensities  [2,3], 
Moreover,  the  APTs  are  transferred  from  one  atom  or  group  of  atoms 
of  a molecule  to  another,  provided  they  occur  in  similar  chemical 
environments  [4],  Those  intensity  parameters  may  be  used,  for 
instance,  to  predict  the  computer-simulated  spectrum  of  new  molecules 
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from  intensity  parameters  transferred  from  known  molecules  [5] . 

This  can  be  advantageously  used  to  predict  the  spectra  of  free 
radicals  or  reactive  intermediate  species  trapped  in  a matrix, 
from  the  study  of  their  infrared  absorption  spectra  and  some 
knowledge  of  their  structure  [5] . The  predicted  spectrum  can 
then  be  compared  with  the  experimental  spectrum  for  qualitative 
identification  of  those  unstable  species.  This  can  be  done, 
for  instance,  in  a computerized  spectrometer  for  better  conve- 
nience of  simultaneous  processing  and  comparison  of  data.  Besides 
this  important  qualitative  application  of  the  theory,  the  quan- 
titative aspect  is  also  of  tremendous  importance.  As  Person  [4] 
has  pointed  out  and  demonstrated  [6],  it  is  possible  to  use  the 
ab  initio  quantum  mechanically  calculated  or  transferred  APTs  to 
predict  the  integrated  molar  absorption  coefficient  (intensity) 
of  a fundamental  band,  not  only  of  neutral  molecular  species,  but 
also  of  free  radicals  [7]  or  other  intermediate  reactive  species 
trapped  in  a matrix  and  then  determine  its  concentration  using  this 
predicted  parameter  and  the  measured  integrated  band  intensity. 

This  potential  analytical  application  of  the  theory  could  be  used 
[4] , for  instance,  for  the  infrared  spectroscopic  measurement  of  the 
concentration  of  free  radical  species  which  dominate  the  chemistry 
of  planetary  atmospheres  [8]  . 

One  of  the  most  interesting  techniques  used  for  selectively 
sampling  trace  species  in  the  stratosphere  is  the  matrix  isolation 
technique  [9] . In  this  technique  an  inert  CO^  matrix  is  used  in 
balloon-born  experiments  to  trap  stratospheric  free  radical  species  , 
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and  after  sample  collection,  the  matrix  is  returned  to  the 
laboratory  for  analysis.  Although  no  measurements  have  yet 
been  reported  [9]  using  the  infrared  spectroscopic  methdd  of 
analysis,  proposed  in  reference  [4] , it  is  expected  that  this 
new  technique  for  analytical  chemistry  [5]  will  lead  to  better 
understanding  of  the  chemistry  of  planetary  atmospheres. 

The  APTs  can  also  be  used  as  a probe  to  investigate  the 
electronic  charge  redistribution  in  molecules  as  a result  of 
molecular  vibrations.  As  shall  be  discussed  in  Chapter  3,  this 
important  piece  of  information  should  be  very  helpful  in  establish- 
ing the  proper  conditions  for  transferability  of  APTs. 

We  chose  the  HOX  series  of  molecules  for  many  reasons.  Our 
interest  in  predicting  the  infrared  spectra  of  matrix  isolated  HO^ 
and  HOF  is  to  show  the  importance  of  this  new  spectroscopic  technique 
for  measurement  of  the  concentration  of  unstable  intermediates  and 
free  radical  species  in  inert  matrices.  The  intensity  parameters 
for  species  such  as  HOF,  HO^  and  HOC£  are  not  easily  obtained  in 
the  gas  phase, and  here  we  think  that  ab  initio  calculations  will 
play  an  important  role  as  an  additional  and  complementary  tool  in 
this  new  developing  area  of  research.  We  also  wish  to  point  out 
that  our  studies  show  for  the  first  time  that  ab  initio  quantum 
mechanical  calculations  can  be  used  to  predict  successfully  the 
relative  intensities  of  unstable  molecules  or  of  free  radical  species 
trapped  in  a matrix  and  indicate  the  importance  of  such  predictions 


using  ab  initio  SCF  calculations. 
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Many  studies  have  been  made  on  the  determination  of  the 
experimental  infrared  spectra  of  matrix- isolated  species  as 
reported  in  the  vast  literature  of  this  subject.  Here  we  think 
it  would  be  interesting  to  use  the  experimental  spectra  parameters 
to  simulate  with  a computer  the  predicted  spectra  of  the  matrix- 
isolated  species.  From  comparison  with  the  experimental  relative 
intensities  the  concentration  of  those  species  in  the  matrix  (if 
the  matrix  pathlength  is  measured)  could  be  obtained  as  well  as 
confirmation  of  the  assignment  of  the  bands  using  the  relative 
predicted  intensities  as  an  additional  criterion.  Moreover,  the 
presence  of  absorption  in  the  spectrum  from  dimers,  trimers,  etc. 
would  be  easily  identified.  This  suggests  another  important 
application  of  the  theory  in  studies  of  inter-  and  intramolecular 
hydrogen  bonded  complexes. 

From  the  molecules  studied  here  (matrix-isolated  H00  and  HOF, 
and  gas  phase  HOC$,) , two  are  of  particular  importance  in  the  chemistry 
of  the  earth's  stratosphere:  H00  and  HOC$,.  Those  species  have  been 

the  subject  of  several  spectroscopic  studies  during  the  past  two  years. 
This  renewed  interest  is  due  in  part  to  the  fact  that  both  species  are 
intermediate  in  the  cycle  of  the  chlorine  catalyzed  decomposition 
of  ozone  in  the  earth's  stratosphere  according  to  the  following 
reactions  [9,10]: 


Cl  + 03  = C$,0  + 02 

(la) 

C$,0  + 0 = C$,  + 02 

(lb) 

°3  + 0 - 2°2 

(lc) 

C$,0  + H0_  = H0C$  + 0„ 
2 2 

(2) 

T_  ^ \ 

hoc$_,;==  ho  + cl 

330nm 

(3) 
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Reactions  (la)  and  (lb)  represent  the  cycle  of  chlorine 
catalyzed  decomposition  of  0^ , while  reaction  (2)  represents  a 
possible  synthesis  for  H0C£  in  the  stratosphere.  It  has  been 
pointed  out  [11]  that  reaction  (2)  is  an  important  removal 
process  for  the  CSLO  species. 

The  rate  constant  of  reaction  (2)  has  been  measured  by  two 

-12 

independent  groups  [12,  p.  2926],  and  a value  of  3.8  + 0.5  x 10 
cm^  molecule  1 s 1 was  reported  by  Reimann  and  Kaufman  [12] . 

This  large  rate  constant  implied  that  H0C£  could  be  an  important 
inert  chlorine  reservoir  in  the  stratosphere  if  its  photo- 
dissociation rate  at  wavelengths  longer  than  300  nm  and  the  reactions 
with  0 and  OH  were  slow  processes  [10] : 

HOCS,  + OH  = CS-0  + H20  (3) 

H0C£  + 0 = C$,0  + OH  (4) 

HOC£  + hV  = c£  + OH  (5) 

The  only  possible  breakdown  of  the  hypothesis  that  HOC#-  is  an  inert 
chlorine  reservoir  in  the  stratosphere  occurs  if  this  species 
would  be  rapidly  photolyzed  at  the  wavelength  region  near  or 
beyond  300  nm.  This  region  is  important  because  at  an  altitude 
of  30-40  km,  light  with  wavelengths  shorter  than  300  nm  is  almost 
totally  attenuated  due  to  the  absorption  by  ozone.  There  had  been 
until  recently  [13]  considerable  controversy  about  the  extent  of 
absorption  and  subsequent  photodissociation  by  HOC#,  in  this  critical 
region.  The  theoretical  and  experimental  studies  of  the  UV  photo- 
dissociation of  H0C&  led  to  different  predictions,  as  did  different 
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experimental  studies.  The  theoretical  SCF-CI  studies  of  the  UV 
photodissociation  of  HOC£  was  done  by  Jaffe  and  Langhoff  [11] . 

They  found  that  the  calculated  absorption  spectrum  of  HOC£ 
exhibits  only  a single  peak  at  220  nm  with  no  significant  ab- 
sorption predicted  beyond  330  nm.  This  theoretical  study  was 
consistent  with  experimental  studies  made  by  Hisatsune  and  by 
Timmons  [14] , but  did  not  agree  with  experimental  studies  by 
Fergusson,  Slotin  and  Style  [15] , by  Jaffe  and  DeMoore  [16] , 
by  Molina  and  Molina  [10]  and  by  Knauth,  Alberti  and  Clausen 

[17].  In  Table  1-1  we  present  several  values  for  the  extinction 

-1-1  1 1 1 " 

coefficient  (£)  in  £ mole  cm  for  the  (1)  A -*■  (1)  A tran- 
sition at  320  nm.  (Some  of  the  measurements  were  made  at  slightly 
different  wavelengths.) 

The  discrepancies  among  the  different  experimental  values 
reported  by  Fergusson  and  DeMore  seem  to  indicate  that  other  species 
besides  H0C£,  H20  and  c£20  were  present  to  different  extents  in 
their  reaction  systems.  The  discrepancies  particularly  at  320  nm 
may  be  explained  if  one  assumes  thermal  and/or  photochemical  de- 
composition of  C£2<3  to  different  extents  in  those  experimental 
works,  combined  with  the  absorption  of  c£20  at  294  nm  which 
probably  have  not  always  been  taken  into  account  to  the  same  extent 
in  the  earlier  experiments  [11] . It  is  well  known  that  c£2  also 
absorbs  at  330  nm  [18]  and  with  a line  shape  similar  to  the  320  nm 
band  observed  in  the  HOC£  experiments . Therefore , it 
seems  that  the  experimental  work  of  Fergusson  is  definitely 
wrong  and  that  the  results  of  DeMoore  at  330  nm  are  probably  over- 
estimated by  factor  of  two. 
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1 ' l " 

EXTINCTION  COEFFICIENTS  FOR  THE  (1)  A -*■  (1)  A TRANSITION 
OF  HOC£  (VAPOR)  NEAR  310  nm 


X 

nm 

e 

(£  mole~lan~l) 

Study 

320 

282 

Fergusson,  Slotin  and  Style, 
(experimental) 

ref.  [15] 

320 

36 

Jaffe  and  DeMoore,  ref.  [16] 

( experimental) 

310 

25.4 

Molina  and  Molina,  ref.  [10] 

( experimental) 

310 

15.7 

Molina,  Ishiwata  and  Molina, 
(experimental) 

ref.  [13] 

310 

16.3 

Knauth,  Alberti  and  Clausen, 
( exp  er imental ) 

ref.  [17] 

310 

1.0 

Timmons,  ref.  [14]  (experimental) 

310 

- 

Hisatsune,  ref.  [14]  (experimental) 

310 

- 

Jaffe  and  Langhoff,  ref.  [11 

] ( theoretical ) 
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All  of  this  controversy  about  the  photochemical  decomposition 
of  HOC&  has  apparently  finally  been  settled  by  the  very  recent 
study  of  Molina,  Ishiwata  and  Molina  [13] . They  have  used  a 
tunable  dye  laser  at  310  nm  to  produce  photodecomposition  of  HOC£ 
measured  by  laser- induced  fluorescence.  It  is  clear  from  this 
study  that  HOC&  does  undergo  photodissociation  by  light  at  310  nm, 
and  so  it  cannot  be  a sink  for  CZ  atoms  in  the  stratosphere. 

Nevertheless,  there  have  still  not  been  any  measurements 
of  the  concentration  of  H0C&  in  the  stratosphere,  and  we  expect 
that  our  work  will  provide  an  alternate  analytical  tool  for 
future  infrared  spectroscopic  measurements  of  concentrations  of 
those  species  in  the  earth's  stratosphere. 


CHAPTER  2 

DESCRIPTION  OF  INTENSITY  CALCULATIONS 
Normal  Coordinate  Analysis 

The  calculation  of  infrared  intensities  using  the  atomic 
polar  tensor  (APT)  formalism  originally  introduced  by  Biarge, 
Herranz  and  Morcillo  [1]  requires  an  a priori  solution  of  the 
vibrational  eigenvalue  problem  represented  by 

G F L = L A (1) 

and  therefore  a complete  normal  coordinate  analysis. 

The  mathematical  approach  used  here  to  solve  the  vibrational 
eigenvalue  problem  is  the  well-known  FG  matrix  method  of  Wilson 
[19,20].  This  method  is  well-documented  [21-24],  and  therefore, 
we  shall  outline  only  the  most  prominent  equations  which  are 
related  to  this  method. 

The  formalism  which  we  shall  discuss  is  applicable  to  a 
general  nonlinear  polyatomic  moleaule  containing  N atoms.  For 
a nonlinear  polyatomic  molecule  there  are  3N-6  normal  modes  of 
vibration  (3N-5  if  the  molecule  is  linear) . 

The  first  step  toward  solving  the  vibrational  eigenvalue 
problem  within  the  Bom-Oppenheimer  and  harmonic  approximation 
consists  of  starting  with  a convenient  set  of  coordinates  and 
then  searching  for  a coordinate  transformation  which  allows  the 
molecular  vibrations  to  be  expressed  in  terms  of  3N-6  (or  3N-5) 
independent  normal  harmonic  motions. 
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We  shall  first  form  a 3Nxl  column  vector  which  represents 
the  displacements  in  position  of  the  atoms  in  the  molecule 
with  respect  to  a space-fixed  Cartesian  coordinate  system 


for  N atoms  of  the  molecule. 

This  set  of  Cartesian  displacement  coordinates  is  not 
independent  since  there  are  six  extra  ("redundant")  coordinates 
(3N  instead  of  3N-6)  associated  with  the  rotations  and  with  the 
translations  of  the  center  of  mass  of  the  system.  Therefore  it 
is  convenient  to  choose  a coordinate  system  (for  example,  molecule- 
fixed  internal  displacement  coordinates)  in  which  these  six 
redundant  coordinates  are  removed. 

The  transformation  from  space-fixed  Cartesian  displacement 
coordinates,  X_,  to  molecule-fixed  internal  coordinates,  R,  is 
given  by 

R = B X (3) 

where  B is  a (3N-6)  x 3N  transformation  matrix. 

Taking  advantage  of  any  symmetry  that  the  molecule  may 
have , and  using  the  sets  of  internal  coordinates  as  basis  functions , 
we  apply  the  projection  operator  technique  to  set  up  the  linear 
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combinations  of  the  internal  coordinates  which  form  the  proper 
symmetry  adapted  functions  transforming  according  to  each  of 
the  irreducible  representations  of  the  symmetry  point  group  of 
the  molecule.  This  set  of  symmetry  coordinates  is  orthonormal. 

In  terms  of  these  symmetry  coordinates,  the  G_  and  F 
matrices  of  the  secular  equation  are  block  factored,  where  each 
block  is  associated  only  with  the  subset  of  symmetry  coordinates 
which  transform  according  to  the  same  irreducible  representation 
of  the  molecular  point  group. 

The  relation  between  the  (3N-6)  x 1 internal  coordinate 
column  vector  and  the  (3N-6)  x 1 symmetry  coordinate  column  vector 
is  given  in  matrix  notation  by 

£ = U R (4) 

where  U is  a (3N-6)  x (3N-6)  orthogonal  matrix,  if.  the  redundancies 
have  been  eliminated. 

At  this  point,  it  is  convenient  to  eliminate  any  redundant 
coordiante  [23,25]  by  dropping  it  from  the  new  basis  set,  since 
those  redundancies  will  only  contribute  to  zero  eigenvalues. 

The  G matrix  and  the  symmetrized  G are  given  by 

G = B M_1  B+  (5a) 

G = U G U+  (5b) 

where  M ^ is  a 3N  x 3N  diagonal  matrix  of  the  inverse  masses. 
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given  by 


M 1 = 


' l/ta_ 


1/frt 


1/m 


1/m, 


1/m, 


1/m, 


(6) 


1/m 


N 


Here  G is  a (3N-6)  x(3N-6)  symmetric  matrix  and  is  related  to  the 

vibrational  kinetic  energy  expression,  given  by 

;+  -1  . f 

2T  = R G R = P G P (7) 

G has  units  of  u \ the  overhead  dot  indicates  a time  derivative, 

the  t is  used  here  to  indicate  the  transpose  of  a matrix,  and  P_ 

is  the  momentum  matrix  conjugate  to  the  internal  coordinates,  R. 

Expanding  the  vibrational  potential  energy  function  V(R^) 

as  a power  series  about  the  equilibrium  configuration  and  neglecting 

2 

terms  of  order  higher  than  R (harmonic  approximation),  we  can 


write  the  vibrational  potential  energy  as 

t 


2V  = E E F.,  R.  R.  = R FR 
10  l 3 — 

13 


(8a) 


where 


p = A2.v 
ij  9r.9r. 
i 3 


(8b) 


or  in  terms  of  symmetry  coordinate  as 

F. 


d2  V 


ij  9S.9S. 

i 3 


(8c) 


13 


Here  F, . and  F . . are  the  force  constants,  or  coefficients  in  the 
13  ID 

vibrational  potential  energy  expression  for  the  molecule  in  terms 
of  internal  or  symmetry  coordinates  respectively. 

The  symmetrized  F_  matrix  is  then  given  by 

F = U F U+  (8d) 


where  F is  a symmetric  (3N-6)  (3N-6)  matrix  whose  elements  are  given 

°-l 

by  Eq.  (8b).  The  F matrix  elements  have  units  of  mdynes  A 

The  units  we  have  used  for  the  various  matrices  are  related 
to  the  units  we  have  chosen  for  the  atomic  masses,  the  potential 
energy  and  the  internal  coordinates.  There  are  two  types  of 
internal  coordinates.  Those  involving  bond  stretching  (Ar  J have 

O 

units  of  A.  Those  involving  angle  deformation  (expressed  as  the 

O O 

arc,  r A0)  have  units  of  (rad  x 1 A = A) . In  other  words,  we  have 

O 

weighted  all  the  angles  by  r = 1 A.  Therefore,  the  resulting 

0-1 

force  constants  (see  Eq. (8b))  have  all  units  of  mdynes  A , when 

O 

the  energy  in  that  equation  is  expressed  in  mdynes  A (1  mdyne 

O 

A = 1 aJ  = 1 attojoule).  The  atomic  masses  are  given  in  atomic 
mass  units  (u) . Table  2-1  gives  a summary  of  matrix  size,  notation 
and  units  used  in  this  work.  Notice  in  Table  2-1  that  use  of 
weighted  units  has  the  great  advantage  of  simplifying  the  matrix 
units . 

Two  important  matrix  relations  from  the  normal  coordinate 
analysis  formalism  are 


L L 


t 


G 


(9a) 


and 


L+  F L = A 


(9b) 


TABLE  2-1 

SUMMARY  OF  MATRIX  SIZE,  NOTATION  AND  UNITS  USED  IN  THE 
NORMAL  COORDINATE  ANALYSIS 


Matrix 

Size 

Matrix 

Units 

. a 
Units 

element 

(unweighted) 

(weighted) 

X 

3N  x 1 

x.  , y. , z. 
i i i 

O 

A 

O 

A 

R 

(3N-6) x 1 

R.  (Ar)b 
i 

o 

A 

O 

A 

R. (A0)C 

rad 

o 

A 

O 

o 

S_ 

(3N-6) x 1 

S. (Ar) 
1 

A 

A 

s. (A0) 
3 

rad 

O 

A 

B 

( 3N-6) x 3N 

B. . (Ar) 
13 

- 

- 

B . (A0) 
kl 

°-l 

A rad 

- 

U 

(3N— 6)  x (3N-6) 

U.  . 
13 

0 1d 

F 

(3N-6)x (3N-6) 

F.  . (Ar, Ar) 
n 

mdyne s A 

mdynes  A 

F. . (Ar ,A0) 
13 

mdynes  rad  d 

°-l 

mdynes  A 

F. . (A0,A0) 

jj 

° -2 
mdynes  A rad 

°-l 

mdynes  A 

M_1 

3N  x 3N 

M.  . _1 
11 

-1 

u 

-1 

u 

G 

(3N-6) x (3N-6) 

G.  . (Ar, Ar) 
n 

-1 

u 

-1 

u 

G.  . (Ar,A0) 
13 

-1  s 

u A rad 

-1 

u 

G.  . (A0,A0  ) 
33 

u-1  A"2  rad2 

-1 

u 

L 

(3N-6) x (3N-6) 

L. . (Ar) 

u_1/2 

u-1/2 

I^i  (AS) 

-1/2  °-l 
u A rad 

u_1/2 

G F 

(3N-6) x (3N-6) 

(G  F)  . . 
13 

-1  ^ n"1 

u mdyne  s A 

-1  °-l 

u mdynes  A 

A 

(3N-6) x (3N-6) 

-1  , °-l 

u mdynes  A 

-1  °-l 

u mdynes  A 

Q 

( 3N-  6 ) x 1 

2k 

1/2  ° 
u A 

1/2  ° 
U A 

funits 

used  in  this 

work 

3Ar  indicates  here  a bond  stretching  displacement 


'A0  indicates  here  an  angle  bending  deformation 
1 mdyne  A~^=  lO^dynes  cm“^=  lOONm-^  = lOOJm-2 = 10' 


'18J  A 2 = 1 aJ  A* 
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Upon  substitution  of  Eq.  (8a)  into  Eq.  (8b)  we  can  easily 
obtain  Eq.  (1) 

G F L = L A 


Therefore , we  can  write 

(G  F - A.  1 ) L = 0 (10) 

k — 3N-6  — 

where  is  the  identity  matrix,  A is  the  diagonal  matrix  of 

eigenvalues,  and  L is  a square  matrix  of  the  eigenvectors  of  (G^ F) . 
The  matrices  (G  F)  , L,  I_  and  A are  (3N-6)  x (3N-6)  matrices. 
Equation  (10)  is  equivalent  to  the  linear  equations 


3N-6 

Z [ (G  F)  . . - 6.  . A ]L  = 0 
j=1  13  13  k jk 


(ID 


where  i = 1,  ...,  3N-6  and  6 is  the  Kronecker  delta. 

Equation  (11)  is  a set  of  linear  homogeneous  simultaneous  equa- 
tions which  have  a nontrivial  solution  if  and  only  if 


| (G  F)  . . - 6.  . A,  |=0  (12) 

13  13  k ' 

The  3N-6  roots  (A  ) of  this  secular  equation  (Eq.  (12))  are  the 
eigenvalues  of  (G  F_)  which  upon  substitution  into  Eq.  (11)  gives 
the  column  vectors  of  L,  which  are  the  eigenvectors  of  the  (G  F) 

product  matrix. 

, -1  <»-l 

Here  A^  has  units  of  u mdyne  A and  is  related  to  the 

harmonic  frequencies  (in  cm  "S  by 

Ak  = 10_5Na  4tt  2c20^  = (1302. 776)  ■2ooJ  (13) 

where  is  Avogadro's  number  and  c is  the  velocity  of  light  in  cm  s \ 
Using  Eq.  (5b) , Eq.  (8d)  and  A,  we  can  show  that 

ff  ? t = I A (14) 

where  the  symmetrized  L matrix  is  given  by 


L = U L 


(15) 
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We  now  proceed  to  obtain  the  normal  coordinates,  0,  which 
are  related  to  the  symmetry  coordinates,  S_,  by  the  linear  trans- 
formation 


Q = L"1  S 


(16) 


P-1 


, Since  L_  is  in  block  diagonal  form,  we  can  also  write  for  each 


symmetry  block 


2<s>  = L-1<S)s<s> 


(17) 


(s) 


where  S_  is  a n-dimensional  column  vector  of  symmetry  coordinates, 

which  transform  according  to  the  same  irreducible  representation  s, 

and  span  the  n-dimensional  space  represented  by  the  column  vector 

(s) 

of  normal  coordinates,  £ 

The  inverse  of  the  symmetrized  matrix  L can  be  obtained  from 


Eq.  (9b)  and  Eq.  (15)  and  is  given  by 

L'1  = A_1  r F 


(18) 


The  transformation  from  a column  vector  of  normal  coordinates 
to  the  column  vector  of  space-fixed  Cartesian  coordinates  X_,  is 
given  by 


X = A R = AULQ=/4LQ 


(19) 


where  A is  the  normal  coordinate  transformation  matrix  from  R to  X 
and  is  given  by 

A = M_1  B+  G-1  (20) 

and  A.  is  the  symmetrized  A matrix. 

Using  the  A matrix,  as  given  by  Eq.  (20),  and  also  Eq.  (1)  to 
form  the  (G  ^ L)  product  matrix 

G-1  L = F L A_1  = (L_1)+ 


(21) 
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we  can  then  form  the  A L_  transformation  matrix 

A L = M_1  B+  G_1  L = M_1  B + F L A_1  (22) 

Finally,  the  symmetrized  C4_  L ) product  matrix  is  obtained  from 
Eq.  (15),  Eq.  (21)  and  the  symmetrized  B_  and  A matrices,  and  is 
given  by 

AL  = m'1/?!!'1  (23) 

where 

B_  = U B (24) 

and 

A = A U+  (25) 

-1/2 

The  (A_  L)  product  matrix  has  units  of  u . This  trams  formation 
is  particularly  useful,  as  will  be  seen  in  the  next  section,  when  we 
shall  consider  the  atomic  polar  tensor  formalism  in  different 
molecular  frames. 

Table  2-2  gives  a summary  of  useful  matrix  equations  involved 
in  the  normal  coordinate  analysis  formalism. 

The  algebraic  equations  given  by  Eq.  (11)  and  Eq.  (12)  are 
not  suitable  for  computer  solution  of  the  vibrational  eigenvalue 
problem  of  large  molecules  with  low  symmetry.  Noticing  that  L is 
nonsingular,  we  can  also  formulate  the  vibrational  eigenvalue  problem 
according  to 

L-1  (G  F)  L = A (26) 

that  is,  since  (G[  F)  is  known  and  L_  is  nonsingular,  it  is  then 
possible  to  find  the  matrix  L which  diagonalizes  (G  F_)  . However, 
since  ((3  F_)  is  not  symmetric  and  since  in  general,  the  most  efficient 
methods  for  computer  solution  of  eigenvalue  problem,  such  as  Eq.  (26) 
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TABLE  2-2 

SUMMARY  OF  USEFUL  MATRIX  EQUATIONS  INVOLVED  IN  THE  NORMAL 
COORDINATE  ANALYSIS  FORMALISM 


Using  internal 
coordinates 

Using  symmetry 
coordinates 

R = B X 

S_  = U R 

“1  *f* 

G = B M B ' 

G_=  U G U 

32  V 

lj  3r.3r. 
i 3 

p = 32  v 

ij  3S.3S. 

i 3 

F = U F U+ 

L L+  = G 

L A+  = G_ 

L+  F L = A 

L+  F A = A 

G F L = L A 

G_F_L  = L A 

3N-6 

Z [ (G  F) . . - 6. . A,  ] L.  = Oa  i=l , , 3N-6 
j=1  13  13  k 3k 

3N-6 

ZHG  F)..  - 6 . . A ]L.  = 
j=1  13  13  k 3k 

1 (G  F ) . . - 6 . . A.  1 = 0a 
13  13  k 

| (G  F)  . . - 6.  . A | = i 
' 13  13  k ' 

Q = L-1  R 

g = if 1 s 

L-1  - A-1  L+  F 

L-1  = A-1  F 

X = A L 0 

X = A L Q 

A = M_1  B+  G-1 

A = A U+ 
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TABLE  2-2 — continued 

-1  -1  -1  -1  ^ 
g =flA  l = (L  ) L 


-1  + 

A L = M B FLA 


,-l 


t 

B 


(L  1) 


t 


CT1  = F_L  A 1 L-1=  (A-1)  L-1 
A L = M_1  B+  F L A'1 


-1  t -1 

= M B (L  ) 


L = U L_ 
B = Uf  B 


L = U L 


B = U B 


a 

See  also  alternative  techniques  suitable  for  computer  solution 
of  the  vibrational  eigenvalue  problem  given  by  Schacht Schneider 
in  ref.  [26]  and  by  Overend  in  ref.  [27]. 


20 


involve  diagonalization  of  symmetric  matrices , alternative  methods 
have  been  mostly  used  [26,27] . Among  other  methods,  for  instance, 
one  which  is  used  in  this  work  involves  two  successive  Jacobi 
diagonalizations  of  G and  F.  We  shall  now  briefly  outline  this 
method  [26]. 

The  first  step  is  a simple  Jacobi  diagonalization  of  G (step  I) 
E?  GD  = T_  (27) 

where  D_  is  orthogonal. 

1/2 

Next,  we  compute  D_  £_  and  form  the  symmetrized  fl  matrix 
(step  II) 

W = D r1/2  (28) 

or 

^ F_W=  H (29) 

Equation  (29)  can  also  be  written  as 

r1/2  D+  F D r1/2  = H (30) 


since 


wf  = r1/2  Df 


It  can  be  easily  shown  that  Eq.  (27)  is  equivalent  to 

r"1/2  Df  G D r_1/2  = I_  (31) 

— — — 3N-6 

where  I_,  . is  the  identity  matrix.  After  postmult  ip  lying  both 
— 3N— o 

sides  of  Eq.  (31)  by  H_ , we  obtain 

r_1/2  0+  G_  D r"1/2  B_  = H (32) 

Now  we  substitute  H_,  as  given  by  Eq.  (30),  into  Eq.  (32) 
r ^/2  G D T~1/2  T1/2  F D T1/2  = 


(33) 


21 


Since  £ is  orthogonal,  Eq.  (33)  becomes 

r“1/2  Df  g f d r1/2  = H 


(34) 


Since  H_  is  symmetric,  it  can  be  diagonalized  by  the  Jacobi 
method  (step  III) 

£+££=A  (35) 

Notice  that  after  finding  the  matrix  £ which  diagonalizes  H we 
can  obtain  the  eigenvectors  of  (£  F)  by  simply  observing  that 

£+  r_1/2£+  G_F  D r1/2  C = A (36) 


Therefore,  the  eigenvectors  of  (£  £)  will  be  (step  IV) 

££1/2£=£  or  (37) 

£ £ = £ (38) 

Up  to  this  point,  we  have  considered  only  the  calculation 
of  vibrational  frequencies  and  normal  modes  from  an  assumed  force 
field  F,  the  atomic  masses  and  the  equilibrium  molecular  geometry. 

We  shall  briefly  consider  now  the  reverse  process,  that  is,  the 
calculation  of  the  force  field  and  normal  modes  from  the  observed 
vibrational  frequencies,  the  atomic  masses  and  the  equilibrium 
molecular  geometry  [28] . The  solution  of  this  problem  is  not  a 
straightforward  procedure  [25] . In  a nonlinear  N-atomic 
molecule  there  are  3N-6  normal  modes  of  vibration,  but  (3N-6) (3N-5)/2 
independent  force  constants  in  the  general  harmonic  force  field. 
Obviously  this  number  is  much  larger  them  the  number  of  observed 
frequencies  for  a single  molecule.  For  very  simple  molecules, 
the  general  harmonic  force  field  can  be  obtained  by  using 
experimental  frequencies  from  isotopic  derivatives  or  other  additional 
experimental  data  (such  as  Coriolis  and  centrifugal  distortion 
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constants).  For  most  molecules,  however,  the  number  of  experimental 
data  available  is  usually  smaller  than  the  number  of  force  constants 
to  be  calculated  and  in  order  to  make  the  problem  solvable,  one 
usually  must  assume  that  certain  force  constants  can  be  neglected. 

As  usual  one  starts  with  a simple  model  force  field  (valence 
force  field  or  Urey-Bradley  type)  obtained  by  transfer  of  force 
constants  from  related  molecules  and  makes  use  of  the  overlay 
least  squares  procedure  to  fit  the  observed  frequencies.  Since 
detailed  presentations  of  this  method  have  been  given  elsewhere 
[22,24-26],  we  shall  only  summarize  the  sequential  steps  used  in 
this  iterative  process. 

We  shall  first  start  with  an  initial  force  field  F°,  which  can 
be  obtained,  for  instance,  by  transfer  of  some  of  the  force  constants 
from  related  molecules. 

We  next  solve  the  vibrational  eigenvalue  problem  for  the  eigen- 

I 

values  Ak  s and  eigenvectors  L using  the  initial  force  field,  F°  [28]. 
From  L we  can  compute  the  Jacobian  matrix,  J,  the  elements  of  which 


are  the  partial  derivatives  of  A s with  respect  to  the  force  constants 


It  can  be  shown  [28] , that  the  difference  between  the  observed 
and  calculated  eigenvalues,  6A,  is  related  to  the  first  order  correction 
to  F , according  to 


The  weighted  least  square.s  solution  of  Eq.  (40)  is  then  obtained 


J 


(39) 


6A  = j 6f 


(40) 


(41) 
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where  W is  a diagonal  matrix  of  weighting  factors.  This  first 
order  correction,  <5f  , is  applied  to  the  previous  force  constants 
to  give  a new  matrix  F 

F = F°  + <Sf  (42) 


Since  Eq.  (40)  is  an  approximation  derived  from  the  assumption 
that  6\  is  a linear  function  of  <5f,  the  cycle  must  be  restarted, 
using  now  the  new  matrix  F_  obtained  in  the  previous  stage.  J and 
6_F  are  recalculated  and  the  process  is  repeated  for  each  cycle 
until  the  calculation  converges  to  an  optimum  set  of  force  constants 
giving  the  best  fit  to  the  observed  frequencies. 

The  elements  of  the  F_  matrix  can  be  expressed  in  terms  of  a 
chosen  set  of  force  constants  (j)  by  the  relation  [26] 


l 2 
l 


jk  a 


(43) 


Equation  (43)  cam  be  represented  in  matrix  notation  as 


F = Z £ (44) 

where  Z is  a retangular  matrix  of  dimension  nxm  where  n is  the 

number  of  force  constants  in  the  F matrix  and  m in  the  number  of 

parameters  in  <j).  Upon  diff  erenti'ation  of  Eq.  (44),  we  obtain 

6_F  = Z_  <5_4>  (45) 

Substitution  of  Eq.  (45)  into  Eq.  (40)  gives 

6A  = Z_  6j)  (46) 

hence,  the  weighted  least  squares  solution  of  Eq.  (46)  give  the 

normal  equations 
"t*  *t* 

Z J W 6\  = (Z  J W J z){j  (47) 

Finally,  the  solution  of  Eq.  (47),  <5<j),  is  obtained  provided 


(48) 


that  (Z_  J W J Z)  is  nonsingular 

<54>  = (z+  j1  wjz)_1(z'  jf  w 5A) 

Failures  to  converge  are  in  general  due  to  a poor  model  force 
field,  F° , or  to  singularity  or  near-singularity  in  the  normal 
equations  [26] . 

It  is  important  to  point  out  that  the  number  of  force  constants 
to  be  refined  should  be  much  less  than  the  number  of  observed 
frequencies,  if  meaningful  results  are  to  be  achieved.  Zerbi  [25] 
has  suggested  that  in  order  to  assure  convergence  toward  a 
statistically  well  defined  set  of  parameters,  one  should  do  the 
calculations  only  when  the  ratio  of  number  of  experimental  data  to 
the  number  of  parameters  approximates  4 or  3. 

Intensity  Relations 

The  theory  of  atomic  polar  tensors  (APTs)  was  originally  introduced 

by  Biarge,  Herranz  and  Morcillo  [1].  Since,  then,  it  has  been 

reformulated  [2]  and  used  extensively  by  Person  and  coworkers  for 

the  analysis  of  infrared  intensities  of  fluorine  containing  molecules, 

such  as,  CFx  H4_x  [29],  SF&  and  UFg  [30],  XF  [31,32]  and  the  first 

organic  molecules  of  the  family  of  alkanes  [33],  alkyl  halides  [34] 

and  ketones,  alcohols,  aldehydes  and  ethers  [35]. 

• (a) 

In  this  theory,  a tensor,  P^  , associated  with  atom  a,  gives 
the  contribution  to  the  change  in  the  total  molecular  dipole  moment, 

Zip,  due  to  a small  displacement,  r^/of  the  nucleus  a [1].  It  is 
assumed  that  a linear  dependence  between  Ap  and  r^  exists  for  small 
nuclear  displacements  and  that  the  simultaneous  displacements  of  the 
nuclei  have  an  additive  effect  on  the  total  molecular  dipole  moment 
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—r 

change,  Ap.  Therefore,  we  can  write 

Ap  = E P (a)? 
a=l~r  a 


(49) 


for  N atonic  nuclei  of  the  molecule. 

In  the  space-fixed  Cartesian  coordinate  system,  the  atomic 
polar  tensor  for  the  otth  atom  of  the  molecule  is  given  by 


P 

— x 


(a) 


\ 


3p  /9x  9p  /9y  9p  /3z 

x a x x a 

9p  /9x  9p  /9y  9p  /9z 

y a y J a ^y  a 

9p  /9x  9p  /9y  9p  /9z 

z a z J a a 


(50) 


here  p^,  p , and  p^  are  the  components  of  the  total  molecular  dipole 
moment,  p,  an  the  space-fixed  Cartesian  coordinate  system. 

The  total  molecular  polar  tensor  is  formed  by  juxtaposition  of 


the  individual  atomic  polar  tensors,  P 

— x 


(a) 


P = {P(D  ! .(2) 


, according  to 


1 .(a)  j ! (N) 


—x 


—x 


! PW>  (51) 


—x 


To  transform  the  total  molecular  polar  tensor  P in  space-fixed 

—x 

Cartesian  coordinate  to  normal  coordinate  space  [2] , we  shall  use 
the  (4_ L)  transformation  matrix,  obtained  from  the  normal  coordinate 
analysis,  as  given  previously  in  Eq.  (23) 


= P A L 
-Q  -x 

-1  + 

= P M B F L A 

—x  — 


(52) 


(53) 


here  P^  has  units  of  e,  A is  unitless  and  L_  has  units  of  u (see 

also  Table  2-1). 

The  integrated  molar  absorption  coefficien,  A^,  associated  with 
the  kth  normal  mode  is  related  to  the  square  of  the  absolute  value 
of  the  total  dipole  moment  derivative  with  respect  to  the  kth 
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normal  coordinate  by  [36] 

\ = (NTrdk/3c2)|9j/9Q  |2 


(54) 


here  N is  Avogadro' s number,  d^  is  the  degeneracy  of  the  kth 
fundamental  mode  and  c is  the  velocity  of  light.  When  is 
expressed  in  km/foole  [2] , Eq.  (54)  becomes 
Ak  = 976.6|3p/9Qk|2 

'Spx  \2  + /9e> 


= 976.6 


,9°k 


+ /9pz 

V 3Q,. 


2 


(55) 

(56) 


3p 

where  x 
9Qk 


9p 

ZQ, 


3p 


and  z are  the  components  of  the  dipole 

3Q. 


moment  vector  (3p/3Qk)  along  each  of  the  space-fixed  Cartesian 

-1/2 

axes,  and  each  has  units  of  e u . 

The  3 x (3N-6)  matrix  of  the  dipole  moment  derivatives 
with  respect  to  the  normal  coordinates  are  formed  by  juxtaposition 
of  each  individual  column  vector,  [2] 


P = { p ip  i 

-Q  1 


X1 


-2  ! 


P 

-Q 


(57) 


3N-6 


(58) 


3V32i  •••  3V3Q3S-6\ 

3py/9Sl  ...  3py/3S3Ht6 

3p/38l  •••  3pzi/3e3N-6  I 
These  dipole  moment  derivatives  with  respect  to  the  normal 
or  Cartesian  coordinates  may  be  obtained  by  two  methods.  The  first 
method  involves  the  measurement  of  the  band  intensities  for  each 
fundamental  mode  of  the  molecule  and  requires  an  a priori  knowledge 
of  the  direction  and  sign  of  the  dipole  derivatives,  as  is  seen  by 


inspection  of  Eq.  (54).  The  experimental  P^  matrix  is  then  formed 
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according  to  Eq.  (57)  and  transformed  to  a space-fixed  Cartesi 
coordinate  to  give  the  experimental  polar  tensor  matrix  P . 


This  transformation  is  given  by  [2] 


— x 


an 


P 

— x 


B + P 3 
- -p- 


(59) 


Here  the _Pp3_ product  matrix  is  the  rotation  correction  term 
necessary  to  transform  from  the  molecule-fixed  normal  coordinate 
space  to  the  space-fixed  of  Cartesian  coordinates  [2]. 

In  the  second  method,  which  is  the  main  subject  of  the  present 
work,  the  dipole  moment  derivatives  with  respect  to  a space-fixed 


Cartesian  coordinate  system  are  obtained  by  quantum  mechanical 


calculations.  Once  the  theoretical  APTs  have  been  calculated, 
we  then  form  the  theoretical  P^  matrix  according  to  Eq.  (51) . 
This  matrix  is  then  transformed  from  a space-fixed  Cartesian 
coordinate  to  a normal  coordinate  space  according  to  Eq.  (52) , 
and  the  theoretical  intensities  are  therefore  obtained  by  using 
Eq.  (55). 


Description  of  Quantum  Mechanical  Calculations 
We  have  seen  in  the  previous  section  that  the  problem  of 
calculating  theoretical  intensities  reduces  to  the  calculation  of 
the  dipole  manent  derivatives  with  respect  to  space-fixed  Cartesian 
coordinates.  We  have  used  the  ab  initio  quantum  mechanical  SCF 
method  using  a 4-31G  basis  set  to  obtain  the  APTs  and  infrared 
intensities  for  a series  of  triatomic  molecules  (HOF,  HOO  and  H0C1) . 
Our  calculations  were  performed  using  the  program  GAUSSIAN  70  [37] 
developed  by  Pople  and  coworkers  [38-43],  The  4-31G  basis  set, 
which  has  been  extensively  used  by  Pople  and  coworkers  [40-42] , 
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is  an  extended  basis  of  atomic  orbitals  expressed  as  fixed  linear 
combinations  of  Gaussian  functions.  In  this  basis,  each  inner  atomic 
orbital  is  represented  as  a linear  combination  of  four  Gaussian  functions 
and  each  valence  orbital  is  split  into  an  inner  part  represented  by  three 
Gaussian  functions  and  an  outer  part  described  by  one  Gaussian  function. 
The  contraction  coefficients  and  exponents  of  the  Gaussian  functions  used 
in  the  program  GAUSSIAN  70,  and  therefore  in  this  work,  were  obtained 
by  minimizing  the  energy  calculated  for  the  atomic  ground  state  of  the 
first  and  second  row  atoms  as  reported  in  references  [40-42] . The  SCF 
calculations  are  based  on  the  restricted  Hartree-Fock  formalism  for 
closed-shell  molecules  and  spin-unrestricted  Hartree-Fock  formalism 
for  open-shell  species  [37] . 

To  allow  more  flexibility  of  the  Gaussian  atomic  orbitals  in 
representing  changes  in  orbital  size  in  different  molecular 
environments  these  atomic  orbitals  are  replaced  by  rescaled  forms 
[40] , and  the  £-scale  factors  are  varied  to  minimize  the  total 
molecular  energy.  Only  the  inner  and  outer  parts  of  the  valence 
shell  were  optimized  [40-42] . We  have  used  the  standard  valence- 
shell  scale  factors  taken  as  the  average  value  of  the  optimum 
scale  factors  obtained  from  the  series  of  molecules  studied  [40,41]. 

An  optional  subroutine  has  been  also  provided  in  program  GAUSSIAN  70 
to  optimize  these  scale  factors  for  any  molecular  calculation, 
in  general,  involving  first  and  second  row  atoms.  However,  the 
study  of  Pople  and  coworkers  [40-42]  based  on  a series  of  nine  molecules 
also  indicates  that  the  range  of  optimum  valence-shell  scale  factors 
is  quite  small  which  justifies  the  use  of  standard  scale  factors. 

Also  no  polarization  functions  have  been  introduced  in  our  calculations. 
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A minimum  basis  set,  such  as  STO  nG (n  = 1,  6)  is  also 

available  in  program  GAUSSIAN  70.  However,  it  has  been  shown  [44] 
that  the  STO  nG  sets  are  very  poor  in  predicting  vibrational 
properties,  and  therefore  no  attempts  have  been  made  here  to  use 
the  STO  nG  sets  for  calculation  of  dipole  moment  derivatives.  It 
has  also  been  shown  [44]  that  only  calculations  with  extended  basis 
sets  predict  intensities  within  the  experimental  uncertainties. 

We  have  used  in  our  ab  initio  calculations  the  experimental 
equilibrium  geometry  instead  of  the  theoretical  calculated  geometry. 
Pulay  [45]  has  pointed  out  that  this  should  be  the  recommended 
procedure.  He  has  shown  that  calculation  of  force  constants  at 
the  experimental  geometry  not  only  yields  better  results,  but  also 
requires  less  computational  work.  Notice  also  that  we  are  interested 
in  differences  between  the  dipole  moment  calculated  at  a certain 
nuclear  displacement  and  at  the  equilibrium  configuration,  and 
therefore,  the  use  of  experimental  geometries  is  not  expected  to 
produce  any  substantial  effect  in  the  calculated  dipole  moment 
derivative-s.  Schwendeman  [46]  has  also  pointed  out  that  the 
experimental  equilibrium  geometry  should  be  used  for  quantum  mechanical 
calculation  of  spectroscopic  parameters,  since  use  of  the  theoretical 
calculated  equilibrium  geometry  may  lead  to  errors  in  the  calculated 
spectroscopic  parameters. 

Let  us  now  consider  how  we  calculated  the  ab  initio  APT  elements 

(a ) 

of  the  P^  matrix  given  by  Eq.  (50)  for  the  ath  atom.  The  APT 
(ct) 

elements  of  the  matrix  were  obtained  in  separate  calculations 

by  displacing  the  ath  atom  from  its  equilibrium  configuration  by 
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o o 

bY  +0.02  A and  -0.02  A along  the  x,  y and  z coordinate  axes, 
and  taking  the  average  value  of  the  dipole  moment  derivatives 
calculated  at  plus  and  minus  displacements  for  each  element 
of  the  tensor.  Let  us  consider,  for  instance,  the  calculation 
of  the  dipole  moment  derivatives  for  the  elements  of  the  first 
column  of  the  P^  matrix.  We  have  approximated  the  dipole 


derivatives  as 


3p.  Ap . 

i l 


3x 


a 


Ax 


a 


<y  _ ply  - <y  |p | w > 
_el  ^ el  i el  ^ 1 el  i „ 

Ax 

a 


(60a) 


Here  <V ^ |p  1^^  represent  the  ith  component  (i  = x,  y,  or  z) 
of  the  total  electric  dipole  moment  of  the  molecule  when  the 
cxth  atom  is  displaced  from  its  equilibrium  configuration  by 

O O /v. 

+0.02  A or  -0.02  A,  <¥  |p|¥  >.  represent  the  i (i  = x ,y  ,or  z ) 

ej-  ei  i o o o o o 

component  of  the  total  electric  dipole  moment  of  the  molecule, 

calculated  at  equilibrium  molecular  configuration,  ¥ is  the 

el 

electronic  wave  function  for  the  molecular  ground  state,  p is  the 

t O 

electric  dipole  moment  operator  and  Ax  = x.  - x = +0.02  A 

a i o — 

• • (OL ) 

Similarly , the  second  and  third  columns  of  P were  obtained 

—x 


8P1  .Api  ~ 

Ay„ 


3y  3y 
a 


(60b) 


a 


and 


3p.  Ap.  <y  |p|y  > - <y  |p|y  >. 

1 ~ 1 _ el 1 1 el  i el 1 ^ 1 el  i, 

Az 


3z  Az 
a a 


a 


where  Ay  = y.  - y = +0.02  A and  Az  = z.  - z 
o.  i o — a i o 


(60c) 

O 

+0.02  A. 
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We  may  therefore  represent  the  dipole  moment  derivatives  as 


(61) 


where  j = (x,  y,or  z) . 


Transferred  APTs 


We  shall  consider  now  the  transferability  of  experimental 
APTs  and  how  we  also  predicted  the  infrared  intensities  for  some 
of  the  molecules  studied  here. 

We  have  pointed  out  previously  that  the  APTs  in  normal  coordinate 
space  for  each  atom  of  a molecule  can  be  determined  experimentally. 

We  have  also  indicated  how  the  experimental  P^  matrix  is  trans- 
formed to  the  P^  matrix.  This  transformation  yields  the  intensity 
parameters  (the  experimental  APTs  in  space-fixed  Cartesian  coordinate) 
which  can  be  transferred  from  one  molecule  to  another,  provided 
they  are  transferred  from  chemically  similar  environments.  We  have 
properly  transferred  some  of  the  experimental  polar  tensors  from 
atoms  of  simple  model  molecules,  such  as,  CH^Cl,  ClCN,  CH3°Hf 

etc.,  to  the  terminal  atoms  of  the  molecules  studied  here  (HOF,  H02, 
and  H0C1)  and  calculated  the  polar  tensor  of  the  central  atom 
using  the  null  matrix  relation  [2] 


for  N atoms  of  the  molecule.  (Here  N is  equal  to  3) . The  polar  tensor 
for  the  central  atom  is  not  expected  to  be  transferable  due  to  the 
uniqueness  of  the  chemical  environment  of  the  central  atom  of  a 


(62) 


molecule. 
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The  experimental  APT  of  a particular  atom  of  a molecule 

(which  was  called  here  the  model  molecule)  was  transferred  in  a 

bond  coordinate  system  for  both  atoms  involved  (the  atom  of  the 

model  molecule  and  the  atom  to  which  the  APT  was  transferred) . 

The  bond  coordinate  system  containing  that  particular  atom  was 

then  rotated  to  the  space-fixed  Cartesian  coordinate  system  of  the 

particular  molecule  studied  .[29].  The  APT  for  both  terminal  atoms 

were  then  transferred  in  the  bond  coordinate  system  and  the  polar 

tensor  of  the  central  atom  was  obtained  using  Eq.  (62).  The  P 

—x 

matrix  was  then  formed  according  to  Eq.  (51)  and  transformed 
to  P^  using  Eq.  (52)  as  indicated  previously.  Equation  (55) 
was  then  used  to  calculate  the  intensities  for  each  fundamental 
band  of  the  molecules  studied. 

The  infrared  intensities  as  calculated  from  ab  initio  APTs 
shall  be  compared  with  the  intensities  predicted  via  transferred 
APTs  as  discussed  in  this  section  and  shall  be  the  subject  of  the 
next  chapter. 


CHAPTER  3 

TRIATOMIC  HOX  MOLECULES 
Hydrogen  Hypofluorite 

We  may  now  apply  the  equations  presented  in  Chapter  2 to 

calculate  the  theoretical  intensities  for  HOF.  By  inspection 

of  Eq.  (53)  given  in  Chapter  2,  we  see  that  the  matrices  P , 

— x 

M i B , F , L and  A have  first  to  be  evaluated  in  order  to 
obtain  the  P^  matrix.  Once  the  P^  matrix  has  been  evaluated, 
the  theoretical  intensity  for  each  fundamental  band  of  HOF  is 
then  calculated  using  Eq.  (56) . 

The  normal  coordinate  analysis  was  carried  out  using  the 
programs  WMAT  and  CHARLY  [27]  from  the  Molecular  Spectroscopy 
Laboratory  at  the  University  of  Minnesota.  The  programs  were 
originally  written  for  the  CDC  6400  computer  and  were  revised 
by  Sanchez  [47]  for  use  on  the  IBM  370/165  computer  at  the 
University  of  Florida. 

The  4-31G  calculations  of  dipole  moment  derivatives  were 
carried  out  using  the  QCPE  program  GAUSSIAN  70  [37] . The  intensity 
calculations  were  performed  using  program  PVDTEN,  developed  by 
Newton  [48] . 

Normal  Coordinate  Analysis 

To  calculate  M , B_  , F_,  L_  and  A \ a normal  coordinate  analysis 
was  performed  using  the  equilibrium  molecular  geometry  from  Kim, 
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Pearson  and  Appelman  [49] . Figure  3-1  defines  the  coordinates 
for  HOF . The  structural  parameters  obtained  by  those  authors 

o o 

were  r(OH)  = 0.964  A,  r(OF)  = 1.442  A and  0 (HOF)  = 97.2°  [49]. 

These  parameters  were  obtained  from  the  analysis  of  the  micro- 
wave  spectrum  of  HOF . Using  these  geometrical  parameters  we 
calculate  the  equilibrium  position  vector  for  the  HOF  molecule, 
given  in  Table  3-1.  These  equilibrium  position  vectors  for  all 
of  the  atoms  of  HOF,  were  used  to  form  the  column  vector,  X, 
given  by  Eq.  (2).  This  column  vector,  the  masses  of  the  atoms 
(H  = 1.007825,  O = 15.994910,  F = 18.998405  in  amu) , the  descrip- 
tion of  the  internal  coordinates  and  the  U matrix  were  used  as 
input  to  the  computer  program  WMAT  [27]  to  obtain  the  B and  G 
matrices  (for  a triatomic  molecule,  belonging  to  the  C symmetry 
point  group,  such  as  HOF,  the  U matrix  is  the  identity  matrix, 
the  internal  coordinates  are  therefore  the  symmetry 
coordinates  of  the  molecule).  In  Table  3-2,  we  present  the  B 
matrix  for  HOF . Program  WMAT  also  diagonalizes  the  G matrix 
and  forms  the  W matrix  according  to  Eq.  (28).  The  F matrix  was 
taken  from  Ogilvie  [50].  It  was  obtained  by  fitting  the  observed 
infrared  frequencies  for  HOF  and  DOF,  using  a multiple  regression 
technique  [50].  The  F matrix  is  given  in  Table  3-3.  The  W matrix 
obtained  from  program  WMAT  and  the  F_  matrix  were  then  used  as  input 
to  program  CHARLY  [27]  to  obtain  the  eigenvectors,  L,  and  eigen- 
values, A,  of  the  G F product  matrix.  The  L matrix  and  the  calculated 
vibration  wavenumbers  of  HOF  are  shown  in  Tables  3-4  and  3-5, 
respectively.  Table  3-5  also  includes  the  experimental  values  for 
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Z 


h3 


3-1.  Coordinate  system  and  atom  numbering  for  HOF. 
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TABLE  3-1 

EQUILIBRIUM  POSITION  VECTORS  FOR  HOF (UNITS  ARE  A) 


Atom 

XT  3 

No. 

X 

y 

z 

1 

0.0 

0.0 

0.0 

2 

o 

• 

o 

-1.442000 

0.0 

3 

0.0 

0.120821 

-0.956399 

a 

The 

atom 

numbering  and 

coordinate  system  are 

given  in  Fig.  3-1. 

THE  B MATRIX  FOR  HOF 
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Note 


THE 

TABLE  3-3 

EXPERIMENTAL  F MATRIX  FOR  HOF (UNITS 

ARE  MDYNE 

3276 

0.4999 

0.0 

4999 

0.9330 

0.1000 

0 

0.1000 

7.2055 

-1 


) 


• F = F 
’ 11  r ( OF ) ' 

F = F 
33  r (OH) . 


F12  Fr (OF) , 0 ' F22  F0;  F23  Fr(OH),0  and 

This  F_  matrix  was  taken  from  Ogilvie  [50] . 


Notice  that  the  following  force  constants  given  by  Ogilvie 
are  scaled  by  the  OF  and  OH  bond  lengths  of  HOF,  according 
to:  Fg  = FQ/ri  r2 , F^Q  = and  F^Q  = F^/r^  Where 

r^  and  r^  are  the  OH  and  OF  bond  lengths  respectively  and 

the  prime  indicates  here  the  scaled  force  constants  given  by 
Ogilvie.  In  addition,  the  force  constants  reported  by  Ogilvie 
are  in  N m and  we  have  transformed  those  force  constants  to 

our  units,  using  the  relation:  F.  . (N  m-1)  x 10-2  = F (mdyne T "S . 

13  l] 

Finally,  notice  also  that  P = F in  this  case,  since  the  U 
matrix  is  the  identity  matrix . 
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i 

Note : 


TABLE  3-4 

THE  L MATRIX  FOR  HOF (UNITS  ARE  u_1'2) 

0.009048  0.004638 

0.030387  -1.081596 

1.026951  0.010244 

Here  L_  = L , since  the  U matrix  is  the  identity  matrix 
in  this  case. 
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TABLE  3-5 

OBSERVED  AND  CALCULATED  VIBRATIONAL  WAVENUMBERS 
FOR  HOF (UNITS  ARE  cm-1) 


a 

Mode 

b 

Observed 

0 

Calculated 

_ d 

Percent  error 

v (OH) 

3578.5 

3591.23 

0.36 

V2(A0) 

1354.8 

1357.47 

0.20 

V3 (OF) 

889.0 

890.63 

0.18 

aHere  (OH)  refers  to  the  OH  stretching  mode,  V (A0)  to  the  bending 
bmode  and  V3(OF)  to  the  OF  stretching  band.  2 

The  experimental  vibration  wavenumbers  of  HOF  in  gas  phase  are  given 

in  ref.  [51]  . 
c 

The  wavenumbers  were  calculated  using  the  equilibrium  molecular 
dgecmetry  given  in  ref.  [49]  and  the  force  field  obtained  from  ref.  [50]. 

[V (calc . ) -V ( obs . ) ] X 100/V (obs . ) 
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the  vibrational  wavenumbers  of  HOF  in  gas  phase  as  given  by 
Appelman  and  Kim  [51]. 

Ab  Initio  4-31G  APTs  and  Calculated  Intensities  for  HOF 

We  have  carried  out  an  ab  initio  quantum  mechanical  cal- 
culation, with  a 4-31G  basis,  of  the  atomic  polar  tensors  (APTs) 
for  HOF.  The  APTs  for  the  oxygen,  fluorine  and  hydrogen  atoms 
were  calculated  according  to  the  procedure  described  in  Chapter  2. 
Once  the  polar  tensors  for  all  of  the  atoms  have  been  calculated, 
we  constructed  the  matrix  using  Eq.  (51) . Notice  also,  that 
the  APTs  for  all  of  the  atoms  of  HOF  were  calculated  relative  to 

the  coordinate  system  shown  in  Fig.  3—1.  In  constructing  the  P 

-x 

matrix,  the  juxtaposition  of  each  individual  APT  follows  the 
same  order  of  atom  numbering  shown  also  in  Fig.  3-1.  The  ab  initio 
APTs  for  HOF  are  presented  in  Table  3-6.  The  intensities  for  each 
fundamental  band  of  HOF  calculated  from  these  ab  initio  4-31G  APTs 
are  presented  in  Table  3-7.  The  comparison  with  experimental 
relative  intensities  for  HOF  and  the  simulated  spectrum  will  be 
discussed  in  the  next  section. 

We  have  also  calculated  the  charge,  charge— flux  and  overlap 
contributions  to  the  total  APT  for  HOF,  according  to  the  CCFO  model 
developed  by  King  [52]  and  extended  by  Zilles  [53].  The  charge, 
c^iar?e—  and  overlap  (CCFO)  model  has  been  discussed  in  detail 

in  ref erence[53]  and  here  we  need  only  to  point  out  the  tensor  contribu- 
tion terms  to  the  ab  initio  APT  and  its  physical  significance.  In 
this  model  the  tensor  contributions  result  from  the  exact  quantum 
mechanical  definition  of  the  APT,  within  the  LCAO-MO  formalism. 
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TABLE  3-6 

AB  INITIO  APTs  FOR  HOF (UNITS  ARE  e,  WHERE  1 e = L602  x 10_19C) 


P° 

—x 


PF 

-x 


-0.129 

0 

0 


PH 

—x 


0.497 

0 

0 


0 

0.174 

-0.020 

0 

-0.437 

0.046 

0 

0.263 

-0.026 


0 

0.228 

-0.068 


\ 

/ 


Note:  We  have  used  a 4-31G  basis  set  for  the  ab  initio  calculation 
of  the  APTs  for  HOF . The  P°  tensor  was  calculated  using  the 
null  condition.  The  coordinate  system  is  shown  in  Fig.  3-1. 
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TABLE  3-7 

ABSOLUTE  INTENSITIES  AND  EXPERIMENTAL  RELATIVE 
INTENSITIES  OF  HOF (UNITS  ARE  km  mole-1) 


Band 

V (cm  "S 

A (4-31G) 3 

A (Expl.) b 

V (OH) 

3578.5 

66.1 

(66.1) 

v2 (A0) 

1354.8 

53.6 

39.6 

v3(of) 

889.0 

11.3 

8.7 

Calculated  intensities  from  the  ab  initio  4-31G  APTs  given  in 

.Table  3-6. 
b 

The  experimental  intensity  for  the  OH  stretching  band  was  scaled 
to  fit  the  intensity  of  (OH)  shown  in  parentheses  as  calculated 
from  the  ab  initio  4-31G  APTs.  Hence,  the  intensities  for  the 
bending  and  OF  stretching  modes  are  relative  band  intensities 
(see  discussion  in  the  next  section) . 
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Accordingly,  the  APT  of  a given  atom  A is  defined  as  the 
transpose  of  the  gradient  of  the  dipole  moment  vector  with 
respect  to  the  displacements  of  the  A atom.  As  a result  of 
analyzing  this  gradient,  the  total  APT  is  broken  down  into 
three  tensor  contributions:  the  charge,  charge-flux  and 
overlap  APTs. 

The  charge  tensor  contribution  to  the  total  APT  on  A is 
given  according  to 


here  = Z^  - N^,  where  Q is  the  gross  atomic  charge  on 

atom  A at  equilibrium  molecular  configuration,  Z^  is  the  charge 

on  nucleus  A and  is  the  Mulliken  gross  atomic  population 

on  atom  A [54,55].  Therefore,  the  charge  APT  gives  the  contribution 


to  the  total  APT  from  the  gross  atomic  charge  associated  with  A, 
at  the  equilibrium  molecular  configuration. 

The  charge-flux  tensor  contribution  to  the  total  APT  on  atom 
A is  given  by  [52] 

P^  (CF)  = I R„  ( V*  Qn  )+  a = A,  B,  C,  ...  (64) 


a 


a 


where  R is 
a 


the  position  vector  of  nucleus  a,  Q is  the 

a 


transpose  of  the  gradient  of  the  gross  atomic  charge  on  atom  a and 
the  summation  is  over  all  the  atoms  in  the  molecule.  The  charge-flux 
tensor  reflects  the  electronic  charge  redistribution  resulting 
from  nuclear  displacement.  Expanding  Eq.  (64)  for  a triatomic 
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molecule,  we  obtain  the  charge-flux  contribution  to  the  total 
APT  for  atom  A as  the  following  sum  of  matrix  products 


3Q 

—1 

9y, 


3Q, 


3Q 


B 


3y, 


9y, 


3q 

_2| 

3z_ 


3q 

c 

3z 


(65) 


here  the  derivatives  were  approximated  as  30,  Ao  , where 

l = xj 

3 j A j 

J A J A 

i= (A, B,  or  C)  and  j=(x,y,  or  z) . The  averaged  terms  Ao  were 

l 

calculated  from  the  differences  in  the  Mulliken  gross  atomic  popula- 
tions of  each  atom  from  the  calculation  at  equilibrium  molecular 
configuration,  Q^,  and  those  from  calculations  made  for  plus  and 
Sirius  displacements  of  the  atom  A,  by  0.02  A along  the  x,  y and  z 
axes;  for  example 


aq. 


Aj 


i 1_ 

2 


Qi-gi 

Aj, 


Q.  - Q° 

i i 


Aj. 


(66) 


here  the  + and  - subscripts  indicate  that  those  quantities  were 
calculated  at  + and  - displacements  respectively. 

The  overlap  tensor  (0)  is  a quantum  mechanical  term  which 
has  no  classical  analogue  [52,53]  and  represents  the  contribution 
of  the  hybridization  effects  (among  others)  to  the  changes  in  the 
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dipole  moment  vectors,  as  the  molecule  vibrates.  Here,  this 
tensor  was  calculated  as  the  difference  between  the  total  APT  and 
the  sum  of  the  charge  and  charge-flux  contributions,  according 
to  Eq.  (62),  [53] 

if  (0)  = pj  (TOTAL)  - (pA  (c)  + PA  (CF) ) (67) 

notice  that  similar  expressions  can  be  given  for  all  the  other  atoms 
in  the  molecule. 

The  charge,  charge-flux  and  overlap  contributions  to  the  ab  initio 
4-31G  APTs  for  the  atoms  in  HOF  are  given  in  Tables  3-8,  3-9  and  3-10, 
respectively,  using  the  coordinate  system  shown  in  Fig.  3-1.  We  shall 
discuss  these  results  further  in  Chapter  4. 

Predicted  and  Experimental  Computer-Simulated  Spectra  of  Matrix— 
Isolated  HOF 


The  band  shapes  [56]  in  the  infrared  spectrum  of  the  matrix- 
isolated  HOF  were  represented  by  the  Gaussian  functions 


a(V)—  a exp 
max 


-4(v-Vq)  in  2 


Av 


1/2 


(68) 


Here,  a(V)  is  the  absorbance  at  wavenumber  V,  a is  the  maximum 

max 

absorbance  of  the  band  at  V , the  wavenumber  at  the  band  origin,  and 
^Vl/2  is  the  ful1  width  <in  on  1)  at  half-maximum  (FWHM)  of  the  absorp- 
tion band.  The  maximum  absorbance  is  related  to  the  total  integrated 
intensity  for  this  Gaussian  line  shape  function  by 


/a  (V ) dv  = /it  (4  In  2)  1/2  a Av 

max  1/2 

Band  = 1.0645  a Av  , 

max  1/2 


(69) 
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TABLE  3-8 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APT  ON  THE  HYDROGEN 
ATOM  OF  HOF (UNITS  ARE  e) . 


0 

0.  45 
0 


PH  (CF) 
— x 


s> 


0 

0 

0 


0 

-0.17 

0.01 


0 

-0.08 

0.17 


0 

-0.02 

-0.03 


0 

0.01 

-0.39 


Note:  The  charge  (C) , charge-flux  (CF)  and  overlap  (O)  contributions 
to  the  total  (T)  4-31G  APT  on  the  hydrogen  atom  were  calculated 
according  to  the  CCFO  model  [52,52].  The  coordinate  system  is 
shown  in  Fig.  3-1. 
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TABLE  3-9 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APT  ON  THE  FLUORINE 
ATOM  OF  HOF (UNITS  ARE  e) . 


F 

(c) 


'0.21 

0 

0 


PF  (CF) 
—x 


0 

0 

0 


F 

P (O) 
—x 


0.08 

0 

0 


F 

P (T) 
—x 


•0.13 

0 

0 


0 

•0.21 

0 

0 

0.12 

0.03 

0 

-0.35 

0.01 

0 

-0.44 

0.05 


0 

-0.04 


\ 


0.02  J 


Note:  The  charge  (C) , charge- flux  (CF)  and  overlap  (O)  contributions 
to  the  total  (T)  4-31G  APT  on  the  fluorine  atom  were  calculated 
according  to  the  CCFO  model  [52,53],  The  coordinate  system 
is  shown  in  Fig.  3-1. 
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Note 


TABLE  3-10 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APT  ON 
OXYGEN  ATOM  OF  HOF (UNITS  ARE  e) . 


P°  (C) 
—x. 


-0.25 

0 

0 


0 

-0.25 

0 


0 

0 

-0 


P°  (CF) 
-x 


0 

0 

0 


0 

0.06 

-0.04 


0 

0 

-0 


P°  (0) 
— x 


-0.12 

0 

0 


0 

0.37 

0.02 


0 

0 

0 


P°  (T) 
—x 


-0.36 

0 

0 


0 

0.19 

-0.02 


0 

0 

-0 


: The  charge  (C) , charge- flux  (CF)  and  overlap 

to  the  total  (T)  4-31G  APT  on  the  oxygen  atom 
calculated  according  to  the  CCFO  model [52,53] 
system  is  shown  in  Fig.  3-1. 


THE 


(O)  contributions 
of  HOF  were 
. The  coordinate 
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Equation  (69)  may  also  be  written  as 


/ 


In (I  /I)dV  = 1.0645  a Av,  .. 
o max  1/2 


(70) 


Band 

here  the  absorbance  a(V)  is  defined  as  the  logarithm  to  the  base  e. 


rather  than  log.- (I  /I),  [57] 
10  o 

a (V)  = In  (I  /I) 
o 


(71) 


where  X and  I are  the  intensities  of  the  radiation  before  and  after 
o 

absorption  by  the  sample,  respectively. 

The  integrated  molar  absorption  coefficient  (experimental  absolute 
intensity)  of  an  infrared  absorption  band,  may  be  defined  as  follows [36] 

IdV  (72) 


Vkm  mole-1)  = iooc.Q," " / ln(VI}< 


Band 

Here  A^  is  the  absolute  intensity  of  the  kth  fundamental  band  in 

-1  . -3 

km  mole  , C is  the  concentration  of  the  sample  in  mmole  cm 

(=  mole  l ■*■)  and  H is  the  cell's  pathlength  in  cm.  The  factor  of 

100  was  introduced  in  Eq.  (72)  to  give  A^  in  km  mole  \ 

Equations  (70)  and  (72)  may  be  combined  to  give 


A^(km  mole  "S  = 1.0645  x 10  2 a Av.,  /0/C& 


(73) 


where  A^ , a 


max 


max  1/2 

, C,  and  % have  been  defined  previously. 


Finally,  we  may  now  substitute  a , as  given  by  Eq.  (73) , 

max 


into  Eq.  (68)  to  give 


-1. 


93.94CilA(km  mole  ) 
a(V)  = — exp 


Av 


1/2 


-4 (V-V  ) In  2 
o 


Av 


1/2 


(74) 


The  predicted  and  experimental  computer  simulations  of  the 
matrix  spectra  of  HOF  were  obtained  according  to  Eq.  (74),  except 
that  the  absorbance  values  were  converted  to  transmitances  and 
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plotted  as  a function  of  V.  The  full  widths  at  half -maximum 
(FWHM)  for  the  absorption  bands  were  estimated  from  values 
reported  in  reference [58] as  described  in  Table  3-11.  The 
absolute  intensities  for  the  three  fundamental  modes  of  HOF  were 
calculated  from  the  ab  initio  APTs  of  HOF.  The  relative  intensities 
were  obtained  according  to  Eq.  (73) . The  c£  parameter  was  adjusted 
^til  the  predicted  intensity  of  the  OH  stretching  vibration  at 
3540  cm  fits  the  observed  spectrum  [59],  according  to  Eq.  (73) 

(the  reason  for  choosing  V(0H)  as  the  fitted  band  is  justified  at 
the  end  of  this  section) . The  value  of  this  c£  parameter  is  given 
in  Table  3-11  and  was  used  for  the  computer  simulation  of  both 
spectra  according  to  Eq.  (74).  The  relative  intensities  are  given 
in  Table  3-11  and  the  predicted  and  experimental  spectra  of  the 
matrix- isolated  HOF  are  shown  in  Fig.  3-2.  We  may  now  compare  both 
spectra.  Notice  that  the  experimental  relative  intensities  predicted 
for  the  bending  and  OF  stretching  modes  are  in  remarkable  agreement 
with  the  corresponding  ab  initio  calculated  intensities. 

The  SPECTRUM  program  written  by  Newton  [60]  was  used  to  generate 
the  spectrum  shown  in  Fig.  3-2.  This  program  utilizes  the  Gould 
software  and  plotter  [61]  to  simulate  an  infrared  spectrum  where  each 
band  is  represented  by  a simple  Gaussian  function  and  overlapping 
bands  by  linear  combination  of  individual  Gaussians.  The  input  for 
program  SPECTRUM  includes  the  scale  data  information  of  the  spectrum, 
the  absolute  intensities,  the  position  of  maximum  absorption,  the 
band  limits  and  the  full  width  at  half-maximum  intensity,  for  each 


fundamental  band. 
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TABLE  3-11 

PARAMETERS  USED  IN  THE  COMPUTER  SIMULATION  OF  THE  PREDICTED 
AND  EXPERIMENTAL  INFRARED  SPECTRA  OF  MATRIX- ISOLATED  HOF 


Predicted 

(C&=6 . 86 x 10-4  mmole  cm 


Experimental 
(c£=6.86x  IQ- 4 mmole  cm 


va 

, °~1\ 
(cm  ) 

Av  . b 
1/2 

(cm  i) 

A (abs . )C 
(km  mole  ) 

voa 
(cm  "1) 

Av  b 
1/2 

(cm”i) 

A(rel. ) C 
(km  mole  ) 

d 

a 

max 

3537.1 

6.0 

66.1 

3537.1 

6.0 

(66.1) 

0.71 

1359.0 

4.4 

53.6 

1359.0 

4.4 

39.6 

0.58 

886.0 

4.0 

11.3 

886.0 

4.0 

8.7 

0.14 

Wavenumber  position  of  maximum  absorption  from  the  experimental 

^matrix  spectrum  of  HOF  as  reported  by  Goleb  et  al.  in  ref.  [59]. 

Full  width  at  half-maximum  of  the  band.  Here  we  have  estimated 

the  full  width  at  half-maximum  by  using  the  values  obtained  from 

the  corresponding  infrared  spectrum  of  matrix-isolated  HO  (in  Ar 

matrix)  reported  by  Smith  and  Andrews  in  ref.  [58] . The  infrared 

spectrum  of  HOF  in  am  N2  matrix  is  reported  by  Goleb  et  al.  in 

ref.  [59].  However,  the  FWHM  of  the  absorption  bands  of  HOF  are 

not  reported  in  that  paper. 

.Taken  from  Table  3-7. 
a 

The  maximum  absorbances,  of  the  HOF  bands  as  reported  by  Goleb 

al.  in  ref.  [59],  but  relative  to  base  e.  The  relation  between 
amax  and  \(re1-)  is  given  by  Eq.  (73). 
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The  experimental  observed  infrared  spectrum  of  HOF  isolated 
in  the  nitrogen  matrix  at  8°K  [59]  is  shown  in  Fig.  3-3.  The 
fundamental  bands  of  HOF  are  seen  at  3537.1  cm”1,  1359.0  cm-1  and 
886.0  cm  ^ for  the  OH  stretching,  bending,  and  OF  stretching  modes. 
The  bands  at  3727.2  and  at  1601  cm  1 belong  evidently  to  H20.  The 
small  band  at  3880  cm  ^ is  the  fundamental  of  HF  [59,62],  formed 
from  thermal  decomposition  of  HOF  before  the  codeposition  of  HOF 
with  N2  at  8°K.  The  small  absorptions  around  the  band  at  3537  cm  1 
are  believed  to  be  due  to  (^0)^  polymers  and  H^O  hydrogen 
bonded  to  HOF  or  to  HF  [59] . The  weak  band  at  3637  cm  ^ is  the 
symmetric  stretch  (v^)  of  H2<0.  Finally,  the  small  band  at  1393  cm  ^ 
was  assigned  [59,63]  to  the  bending  mode  of  HOF  hydrogen  bonded  to 
HF.  The  upper  and  lower  spectrum  shown  in  Fig.  3-3  were  recorded 

at  different  concentrations  of  HOF  in  the  N matrix  [59] . 

2 

The  reason  for  choosing  the  V(OH)  band  as  the  one  to  be  fitted 
was  based  upon  the  reliability  of  the  4-31G  basis  set  in  predicting 
the  intensity  of  this  mode.  In  general,  the  agreement  between 
the  predicted  and  experimental  intensities  for  this  mode  is  within 
factor  of  two  or  better,  as  shown  in  Table  3-12  for  a series  of 
molecules  studied  in  our  laboratory. 

Table  3-12  shows  that  the  predicted  intensities  for  CH^OH  and 
HOC1  are  within  a factor  of  1.0  and  1.3  agreement  with  the  corre- 
sponding experimental  intensities.  For  the  and  vibrational 
modes  of  H^O  the  predicted  intensities  are  within  factor  of  2.1  and 
2.3  agreement  with  the  corresponding  experimental  values.  A 
factor  of  two  agreement  with  the  experimental  intensities  is  expected 
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TABLE  3-12 

EXPERIMENTAL  AND  QUANTUM  MECHANICALLY  PREDICTED  INTENSITIES 


FOR  THE  V (OH) 

VIBRATIONAL 

MODE  OF  CH3OH,  H20 

AND  HOC1 

Molecule 

V (OH) 

A (4-31G) 

A (Expl.) 

(gas  phase) 

(cm~l) 

(km  mole--*-) 

(km  mole  ) 

ch3oh 

3667.0 

22. 2a 

22. 0a 

HO 

3633.2 (V, ) 

1.07b 

2.24° 

2 

1 

3725. 7(V  ) 

19. 7b 

44. 6C 

HOC1 

3626.0 

75.3d 

56. 6S 

3Rogers  [35,  p.  118] 
bZilles  [53,  p.  73] 


CQuoted  by  Zilles  [53,  p.  73] 
^This  work,  p.  114 
6 Su  et  al.  [78] 
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from  intensities  predicted  from  ab  initio  SCF  calculations 
using  an  extended  basis  set.  This  level  of  agreement  is  reason- 
able since  it  is  still  not  surprising  to  find  that  a new 
measurement  of  a band  intensity  differs  by  a factor  of  two  from 
earlier  results  reported  in  the  literature. 

Transferred  APTs  and  Predicted  Intensities  for  HOF 

The  intensities  for  HOF  have  been  also  predicted  from 
experimental  APTs  transferred  from  other  molecules.  The  experimental 
fluorine  and  hydrogen  APTs  from  CH^F  [29]  and  H^O  [53]  respectively 
were  transferred  in  a bond  coordinate  system  to  the  molecular 
coordinate  system  of  HOF,  shown  in  Fig.  3-1.  The  oxygen  APT  was 
calculated  using  the  null  condition  given  by  Eq.  (62) . In  the 
first  column  of  Table  3-13,  the  oxygen,  fluorine  and  hydrogen  APTs, 
transferred  from  CH^F  and  H^O  are  shown,  defined  in  the  coordinate 
system  of  Fig.  3-1.  We  have  also  predicted  the  intensities  for 
HOF,  using  the  experimental  APTs  transferred  from  CH^F  [29]  and 
CH^OH  [35].  This  is  shown  in  the  second  column  of  Table  3-13.  For 
comparison,  we  have  also  included  in  Table  3-13,  the  ab  initio 
APTs  for  HOF. 

The  intensities  calculated  from  transferred  and  ab  initio 
APTs  for  HOF  are  shown  in  Table  3-14.  Notice  that  the  intensities 
predicted  either  by  transfer  of  the  APTs  from  CH^F  and  H^O  or  CH^F 
and  CH^OH  are  both  in  poor  agreement  with  the  intensities  calculated 
from  the  ab  initio  APTs  and  also  with  the  relative  experimental 
intensities.  The  transferred  APTs  predict  an  order  of  increasing 
intensities  from  high  wavenumbers  (V^(OH))to  low  (V^(OF)),  precisely 


TABLE  3-13 

TRANSFERRED  AND  AB  INITIO  4-31G  APTs  FOR  HOF  (THE  COORDINATE  SYSTEM  IS  SHOWN  IN  FIG.  3-1  AND  UNITS 
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TABLE  3-14 

CALCULATED  INTENSITIES  FROM  TRANSFERRED  AND  AB  INITIO  APTs 
FOR  HOF (UNITS  ARE  km  mole-1) 


Band 

V (cm  1) 
k 

\a 

\C 

V (OH) 

3578.5 

16.3 

17.0 

66.1 

V2  (A0) 

1354.8 

61.3 

30.0 

53.6 

V3 (OF) 

889.0 

73.5 

81.2 

11.3 

3.  F H 

Calculated  intensities,  using  P and  P from  CH^F  [29]  and  H^O  [53] , 

^respectively  (see  column  1 of  T|fble  3^3) 

Calculated  intensities,  using  J?  and  P_H  from  CH^F  [29]  and  CH^OH  [35]  , 

respectively  (see  column  2 of  Table  3-1)3) 

Calculated  intensities  using  4-31G  APTs  given  in  column  3 of  Table  3-13. 
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the  opposite  of  what  is  observed  experimentally  or  predicted 

quantum  mechanically  (see  Table  3-7  and  Fig.  3-2) . The  predicted 

intensities  for  the  OH  and  OF  stretching  modes  of  HOF,  using 
F H 

P from  CH  F and  P from  HO,  given  in  column  3 of  Table  3-14, 

X J X 4 

are  off  by  factors  of  4 and  7 respectively,  compared  to  the 
corresponding  calculated  intensities  using  ab  initio  APTs. 

Only  the  intensity  predicted  for  the  bending  mode  is  in  reasonable 
agreement  with  the  corresponding  intensity  calculated  from  ab  initio 
APTs.  For  the  intensities  predicted  from  APTs  transferred  from 
CH^F  and  CH^OH,  the  results  appear  to  be  even  worse  (see  column  4 
of  Table  3-14) . These  results  are  not  unexpected,  since  the 
experimental  APTs  in  fluorine  and  hydrogen  transferred  from 
CH^F  and  CH^OH  or  H20  are  substantially  different  from  the  correspond- 
ing ab  initio  APTs  (see  Table  3-13).  We  may  now  compare  the  elements 
of  the  ab  initio  APTs  for  fluorine  and  hydrogen  with  the  correspond- 
ing tensor  elements  transferred  from  CH  F and  H„0  or  CH„OH.  Let 

3 2 3 

us  for  instance,  compare  the  ab  initio  calculated  diagonal  elements 

3p  3p  and  3p  of  the  4-31G  APTs  of  fluorine,  given  in  the 
x y z 

3x  ' 3y  3z 

second  row  of  Table  3-13,  with  the  corresponding  elements  transferred 

from  CH  F and  HO  or  CH  OH.  The  3p  and  3p  tensor  elements  of  the 
3 2 3 x y 

3xF  3yF 

ab  initio  APT  differ  by  a factor  of  two  from  the  corresponding  values 

transferred  from  CH_F.  The  3p  is  off  by  a factor  of  1.5.  It  is 

3 z 

3z 

F 

obvious  that  the  experimental  fluorine  APT  of  CH^F  cannot  be 
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successfully  transferred  to  the  HOF  molecule.  Moreover,  the  3p 

s' 


element  of  the  ab  initio  APT  is  -0.16  e as  compared  to  the 

corresponding  null  element  of  CH^F. 

The  elements  of  the  P tensor  are  also  not  transferable,  as 

—x 

can  be  seen  on  the  third  row  of  Table  3-13.  Here  we  have  to 

consider  the  elements  of  the  P tensor  transferred  from  CH  OH  and 

-x  3 

H^O  separately.  Let  us  consider  first  the  elements  of  the 

hydrogen  APT  of  HO.  The  diagonal  elements  9p  and  3p  differ 
2 x z 


3x. 


3z. 


H H 

by  a factor  of  1.5  and  2.0  respectively,  from  the  corresponding 
elements  of  the  ab  initio  APT.  Only  the  3p  element  of  the 


3y 


H 


transferred  APT  from  H^O  is  in  fairly  good  agreement  with  the 
corresponding  ab  initio  APT. 


For  the  elements  of  the  P tensor  transferred  from  CH  OH, 

-x  3 

the  diagonal  elements  3p  and  3p  differ  from  the  ab  initio  APT 

x z 


3x. 


3z. 


H 


by  factors  of  2.0  and  the  3p  element  differs  by  a factor  of  1.5. 

3yt 


H 


An  important  question  which  may  be  raised  at  this  point  is 

the  following:  Why  was  it  possible  to  transfer  the  experimental 

fluorine  APT  successfully  from  CH^F  to  predict  the  intensities 

of  CF  , CF  H,  CH  F [29]  and  SF.  [64]  but  then  fail  to  predict 
4 3 2 2 6 

the  intensity  of  HOF?  The  answer  is  obviously  due  to  the 
different  nature  of  the  chemical  bonding  of  the  fluorine  atom  in 
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the  original  molecule  compared  to  that  in  the  molecules  to 
which  the  APT  is  being  transferred.  For  the  molecules  studied 
in  reference  [16] , the  electron  distributions  around  the  C-F  bond 
in  CH^F , CF^,  CF^H  and  CH^F^  are  sami^ar'  so  that  the  substitution 
of  the  hydrogen  atoms  in  CH^Fby  fluorine  atoms  does  not  sub- 
stantially change  the  electron  distribution  in  the  C-F  bond. 

For  the  HOF  molecule,  the  situation  is  significantly  different. 

Here,  the  electron  distribution  in  the  X-F  bond  is  different  in 

CH^F  and  HOF  (for  one  reason,  due  to  the  change  in  symmetry)  and 
, F 

the  tensor  may  not  be  expected  to  be  very  transferable. 

It  is  possible  therefore,  to  generalize,  concerning  the 
transferability  of  APTs,  by  stating  that  the  APTs  are  expected  to 
be  transferred  frcm  one  molecule  to  another  providing  the  electron 
distribution  around  the  bond  to  which  the  APT  is  being  trans- 
ferred is  similar  to  that  in  the  original  molecule.  As  a consequence 
of  this  rule,  we  may  also  say  that  within  a series  of  chemically 
similar  molecules  the  APTs  are  expected  to  be  transferable,  since 
the  chemical  environment  near  the  bond  to  which  the  APT  is  being 
transferred  does  not  change  substantially  from  one  molecule  to 
another. 

In  conclusion,  we  should  expect  to  be  able  to  transfer 
successfully,  for  instance,  the  experimental  hydrogen  APT  from 
CH^F  to  CH^F^  or  CHF^  molecules,  but  we  should  not  expect  to  be 
able  to  transfer  the  experimental  hydrogen  APT  frcm  HOF  to  predict 
successfully  the  intensities  of  HOC&  or  HOBr  molecules.  Here 
we  can  expect  that  the  change  from  F to  Ci  to  Br  will  have  a sub- 
stantial effect  in  changing  the  electron  distribution  near  the  OH 


bond. 
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Hydroperoxyl  Radical 

The  hydroperoxyl  radical  HO.,  is  an  important  inter- 
mediate species  in  the  chemistry  of  the  earth's  stratosphere 
[9] . Direct  observation  of  this  species  was  first  accomplished 
by  Foner  and  Hudson  in  1953  [65]  using  mass  spectrometric 
technique  and  by  Milligan  and  Jacox  in  1963  [66]  using  infrared 
spectra  matrix  isolation  technique.  Since  1972,  the  number  of 
spectroscopic  studies  of  HO.,  has  increased  steadily  [58,67-75]. 

The  renewed  interest  in  this  species  was  mainly  due  to  the  fact 
that  it  takes  part  in  several  cyclic  processes  which  remove 
ozone  from  the  upper  and  lower  stratosphere  of  the  earth  [9] . 
Moreover  it  has  long  been  postulated  as  an  important  inter- • 
mediate  in  numerous  hydrogen-oxygen  reactions  (see  for  example, 
references  cited  in  [68] ) . 

Normal  Coordinate  Analysis 

The  molecular  geometry  of  HO,,  in  its  ground  electronic  state 

2 " 

( A ) is  well  established,  both  by  experiment  [71,74]  and  by  ab 
initio  SCF  calculations  [75] . The  structural  parameters  of  HO^ 
are  given  in  Table  3-15,  as  reported  by  different  research  groups. 
The  agreement  between  them  is  remarkably  good.  In  this  work  we 
have  used  the  geometrical  parameters  reported  by  Beers  and  Howard 
[71]  as  given  in  Table  3-15.  From  these  geometrical  parameters, 
we  calculate  the  position  vectors  for  the  HO^  molecule , given  in 
Table  3-16,  using  the  coordinate  system  shown  in  Fig.  3-4.  The 
column  vector  X_  given  by  Eq.  (2)  was  then  formed  using  the  data 
presented  in  Table  3-16.  This  column  vector,  the  masses  of  the 
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STRUCTURAL 

TABLE  3-15 
PARAMETERS  OF  H02 

2 " 
(A  ) 

r (OH) a 

r (00) 3 

0 (HOO) 

References 

0.977 

1.335 

104.1° 

Beers  and  Howard,  ref. 
[71] (microwave) 

0.975 

1.329 

104.0° 

Tuckett,  Freedman  and 
Jones,  ref.  [74] (near 
infrared) 

0.968 

1.384 

106.8° 

Liskow,  Schaefer  III  and 
Bender,  ref.  [75]  (ab 
initio  SCF  calculation) 

0.973 

1.458 

104.6° 

Liskow,  Schaefer  III  and 
Bender,  ref.  [75]  (ab 
initio  SCF-CI  calculation) 

a 

Bond  length 

in  units  of 

O 

A. 
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TABLE  3-16 


EQUILIBRIUM 

POSITION  VECTORS 

FOR  H02 (UNITS  ARE  A) 

a 

Atom  N 

o 

X 

Y 

z 

1(0, central) 

0.0 

0.0 

0.0 

2 (0, outer) 

0.0 

-1.335000 

0.0 

3(H) 

0.0 

0.238012 

-0.947565 

aThe  atom  numbering 

and  coordinate 

system  are  given 

in  Fig.  3-4. 
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Z 


Hj 


Fig.  3-4.  Coordinate  system  and  atom  numbering  for  HO^. 
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atoms  (H  = 1.007825,  0 = 15.994910  in  amu) , the  description 
of  the  internal  coordinates  and  the  U matrix  were  used  as 
input  to  program  WMAT  [27]  to  obtain  the  and  G matrices. 

The  B matrix  for  HC>2  is  presented  in  Table  3-17.  The  G 
matrix  was  then  diagonalized  by  program  WMAT  to  form  the  W 
matrix,  according  to  Eq.  (28) . We  used  the  force  field  reported 
by  Ogilvie  [50] , which  was  obtained  by  fitting  the  observed 
infrared  frequencies  for  eight  isotopic  derivatives  of  the  HC>2 
species,  using  a multiple  regression  technique.  The  F matrix 
is  given  in  Table  3-18.  The  W matrix  obtained  from  program 
WMAT  and  the  F matrix  were  used  as  input  to  program  CHARLY  [27] 
to  obtain  the  eigenvectors,  L,  and  eigenvalues  A,  of  the  G F 
product  matrix.  The  L matrix  and  the  calculated  vibration 
wavenumbers  of  HC>2  are  shown  in  Tables  3-19  and  3-20,  respectively. 

In  Table  3-20  is  also  given  the  experimental  vibration  wave- 
numbers  of  H02 , obtained  from  H02  isolated  in  an  argon  matrix,  as 
reported  by  Ogilvie  in  reference [50] . 

Ab  Initio  4-31G  APTs  and  Calculated  Intensities  for  H02 

The  atomic  polar  tensor  (APTs)  for  H02  were  quantum  mechanically 
calculated  using  a 4-31G  basis  set.  The  details  of  this  calculation 
have  been  previously  discussed  in  Chapter  2.  These  APTs  were 
calculated  with  respect  to  the  coordinate  system  shown  in  Fig.  3-4. 
The  P^  matrix  was  then  constructed  according  to  Eq.  (51) , noticing 
that  the  juxtaposition  of  each  individual  APT  in  this  equation, 
follows  the  same  order  of  atom  numbering  shown  in  Fig.  3-4.  The 
ab  initio  APTs  for  HC>2  are  given  in  Table  3-21.  The  APT  for  the 


TABLE  3-17 
THE  B MATRIX  FOR  HO, 
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TABLE  3-18 

THE  EXPERIMENTAL  FORCE  CONSTANT  MATRIX,  F,  FOR  HO, 
o _2 (UNITS  ARE0MDYNE  &-1) 

(laJ  A =1  mdyne  A”-*-  = 100  N m ) 


5.7186 

0.3870 

0.1000 

0.3870 

0.9870 

0.0500 

0.1000 

0.0500 

6.5817 

Note : 


F11  Fr(Oor  F12  Fr(OO),0'  F13  Fr (00) ,r (OH) 1 

F22  ‘ V F23  * Fr  (OH), 9 ' F33  = Fr  (OH)  ' NOtice  that  4116 

following  force  constants  given  by  Ogilvie  [50]  are  scaled 

by  the  00  and  OH  bond  lengths  of  HO  , according  to:  Fg=  Fg/r  r , 

F 0=  F ft/r  and  F Q=  F Q/r_  , where  r,  and  rn  are  the  OH 
r1t3  r^0  1 r^0  r20  2 1 2 

and  00  bond  lengths  respectively  and  the  prime  indicates  here 
the  unsealed  force  constants  given  by  Ogilvie.  Notice  also 
that^the  force  constants  reported  by  Ogilvie  are  in  units  of 
N m . Here  F = F , since  the  U matrix  is  the  identity  matrix. 
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1 
Note : 


TABLE  3-19 

THE  L MATRIX  FOR  H02 (UNITS  ARE  u /2) 

-0.014835  0.031959 

-0.043829  -1.091023 

1.027013  0.000375 

Here  = L,  since  the  U matrix  is  the  identity  matrix  in 
this  case. 
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TABLE  3-20 

OBSERVED  AND  CALCULATED  VIBRATIONAL  WAVENUMBERS 
FOR  HO  (UNITS  ARE  cm  *) 


Mode3 

Observed*3 

0 

Calculated 

d 

Percent  error 

V (OH) 

3410.50 

3432.58 

0.65 

v2  (A0) 

1389.00 

1399.71 

0.77 

v3 (00) 

1101.00 

1084.51 

-1.50 

Here  (OH) , V (A0)  and  (00)  refer  to  the  OH  stretching,  the  bending 
j_mode  and  the  OO  stretching  band,  respectively. 

The  experimental  vibration  wavenumbers  from  HO  isolated  in  an  argon 
matrix  are  given  in  refs.  [50],  [66]  and  [58].  We  have  used  the  data 

Reported  by  Ogilvie  in  ref.  [50]. 

The  wavenumbers  were  calculated  using  the  equilibrium  molecular  geometry 
of  HO  given  in  ref.  [71]  and  the  force  field  obtained  from  ref.  [50]. 
[V(calc.)  - V (obs . ) ] x 100/V(obs.). 
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TABLE  3-21 

AB  INITIO  APTs  FOR  HO^ (UNITS  ARE  e,  WHERE 
le  = 1.602  x 10-19  C) 


pOU> 

— x 


(-0.433 

: 


o 


-0.697 

0.0 


pO(2) 

-x 


0 

0.393 

0.015 


PH 

—x 


0 

0.304 

-0.015 


Note:  A 4-31G  basis  set  was  used  for  the  ab  initio  calculation  of 
the  APTs  for  HO^  • The  tensor  was  calculated  using 

the  null  condition.  The  coordinate  system  is  shown  in  Fig. 
3-4. 


75 


central  atom  was  calculated  from  the  null  condition  [Eq.  (62)]. 

The  matrix  was  transformed  from  the  space-fixed  cartesian 
coordinate  system  to  the  matrix  in  normal  coordinate  space, 
using  Eq.  (53)  and  the  intensities  for  each  fundamental  band  of 
H02  ware  calculated  using  Eq.  (56) . 

We  have  also  calculated  the  charge,  charge-flux  and  overlap 
contributions  to  the  ab  initio  APTs  in  H02 , according  the  the 
CCFO  model  [52,53]  using  Eq.s  (63),  (64)  and  (67),  respectively. 
These  APTs  are  given  in  Tables  3-22,  3-23  and  3-24.  For  all  of 
these  terms  the  coordinate  system  is  as  shown  in  Fig.  3-4.  The 
charge,  charge-flux  and  overlap  contributions  to  the  APT  for  HO 

2 

shall  be  discussed  in  Chapter  4,  where  we  shall  compare  them 
for  all  of  the  molecules  studied  in  this  work. 

For  the  intensity  calculations  we  have  used  program  PVDTEN 
developed  by  Newton  [48] . The  intensities  for  each  fundamental 
band  of  H02  calculated  from  the  ab  initio  4-31G  APTs  are  given  in 
Table  3-25.  Also  shown  in  that  table  are  the  absolute  intensities 
calculated  by  Komornicki  and  Jaffe  [ 76  ] with  the  GRADSCF  program 
[77]  using  a contracted  Gaussian  basis  of  triple-zeta  accuracy  plus 
double  polarization  (TZPP) . We  have  also  included  in  Table  3-25 
the  experimental  intensities  relative  to  the  intensity  of  the  OH 
fundamental  stretching  mode  V(OH)  of  H02,  obtained  from  the 
experimental  spectrum  of  H02  in  an  Ar  matrix,  reported  by  Smith 
and  Andrews  in  reference  [25] , and  shown  in  Fig.  3-5. 

It  is  of  interest  now  to  compare  our  intensities  calculated 
with  a 4-31G  basis  set  with  those  predicted  by  Komornicki  and  Jaffe 
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Note : 


TABLE  3-22 

CHARGE,  CHARGE-FLUX, OVERLAP  AND  TOTAL  APTs  ON  THE 
HYDROGEN  ATOM  OF  H02 (UNITS  ARE  e) 


PH  (C) 
—x 


0 

0.44 

0 


PH  (CF) 
— x 


0 

■0.12 

-0.02 


PH  (0) 
—x 


PH  (T) 
-x 


0 

-0.02 

0.004 

0 

0.  30 
-0.02 


The  charge  (C) , charge-flux  (CF)  and  overlap  (O)  contributions 
to  the  total  (T)  4-31G  APT  on  the  hydrogen  atom  were  calculated 
according  to  the  CCFO  model  [52,53],  The  coordinate  system  is 
shown  in  Fig.  3-4. 
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TABLE  3-23 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APTs  ON  THE  OUTER 
OXYGEN  (2)  ATOM  OF  HO  (UNITS  ARE  e) 


p0(2)  (C) 

—x 


0 

-0.04 

0 


P0(2)(T) 

-x 


0 

0.39 

0.02 


Note:  The  charge  (C) , charge-flux  (CF)  and  overlap  (O)  contributions 
to  the  total  4-3 1G  APT  on  the  oxygen  (2)  atom  were  calculated 
according  to  the  CCFO  model  [52,53].  The  coordinate  system 
is  shown  in  Fig.  3-4. 
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TABLE  3-24 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APTs 
ON  THE  CENTRAL  OXYGEN  (1)  ATOM  OF  HO 
(UNITS  ARE  e)  2 


0 

-0.23 

-0.003 


0 

-0.03 

0.01 


0 

-0. 66 

0.01 


Note:  The  charge  (C) , charge-flux  (CF)  and  overlap  (O) 

contributions  to  the  total  (T)  4-31G  APT  on  the 
oxygen  (1)  atom  were  calculated  according  to  the  CCFO 
model  [52,53].  The  coordinate  system  is  shown  in 
Fig.  3-4. 
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TABLE  3-25 

ABSOLUTE  INTENSITIES  AND  EXPERIMENTAL  RELATIVE  INTENSITIES 
OF  HO  (UNITS  ARE  km  mole-1) 


Calculated  Relative  experimental 


Band 

V (cm  d) 

A (4-31G) 3 

A (TZPP)b 

c 

A 

Ad 

V (OH) 

3410.5 

47.6 

67.5 

47.6 

67.5 

v2 (A0) 

1389.0 

78.5 

53.2 

148.3 

210.1 

V3 (00) 

1101.0 

25.9 

46.  3 

91.3 

129.4 

a 

Calculated  intensities  using  the  4-31G  APTs  given  in  Table  3-21 
k (this  work) . 

Calculated  by  Komornicki  and  Jaffe  and  reported  in  ref. [76]  using 

a triple-zeta  plus  double  polarization  (TZPP)  basis  set  of  Gaussian 

functions  contracted  to  (593p2d/392p) , for  the  calculation  of  the 

dipole  moment  derivatives, 
c . 

Experimental  intensities  relative  to  A(4-31G).  Relative  experimental 
intensities  scaled  to  fit  the  intensity  of  the  V(OH)  band  calculated 
dby  4-31G  APTs. 

Experimental  intensities  relative  to  A (TZPP).  Relative  experimental 
intensities  scaled  to  fit  the  intensity  of  the  V (OH)  band  calculated 
by  TZPP  basis  set. 
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using  the  TZPP  basis  and  also  compare  both  calculated  intensities 
with  the  relative  experimental  intensities  of  HO^.  The  SCF  cal- 
culations reported  by  Komornicki  and  Jaffe  predict  the  most  intense 
band  to  be  associated  with  the  OH  stretching  mode  which  is  not 
observed  experimentally  [58] . Our  SCF  calculations  predict  the 
bending  mode  to  be  the  strongest  band  of  HO^  in  agreement  with  the 
experimental  relative  intensities.  Our  absolute  intensities  are 
also  consistent  with  the  experimental  observations  of  Jacox  and 
Milligan  [67]  who  pointed  out  that  the  bending  mode  of  HC>2  is 
definitely  stronger  than  the  00  stretching  band.  Let  us  now 
compare  the  absolute  intensities  calculated  with  the  4-31G  and 
TZPP  bases  with  the  corresponding  relative  intensities  from 
experiments.  The  intensity  pattern  for  the  H00  bending  mode 
and  00  stretching  mode  Calculated  with  a 4-31G  basis  and  using 
the  APT  formalism  is  good,  when  compared  with  the  corresponding 
experimental  pattern.  We  have  apparently  underestimated  the  rela- 
tive intensity  only  for  the  00  stretching  mode.  However,  the 
intensity  pattern  calculated  with  the  TZPP  basis  does  not  match 
the  experimental  pattern  so  well  (see  Fig.  3-7  in  the  next  section) 
since  it  predicted  all  bands  to  have  about  the  same  relative  in- 
tensity in  contrast  to  the  variations  found  experimentally. 

Although  the  relative  intensities  from  this  calculation  appear  to 
agree  with  the  experimental  values  within  factors  of  2-4,  we  believe 
that  our  calculated  intensities  are,  in  general,  more  consistent 
with  the  experimental  observation  than  are  those  predicted  with 
the  TZPP  basis  set.  It  is  possible  that  the  TZPP  basis  set  [ 76  ] 

is  not  well-balanced  [55]  leading  to  some  error  in  predicting  dipole 
moment  derivatives,  as  shall  be  discussed  in  more  detail  later. 
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Predicted  and  Experimental  Computer- Simulated  Spectra  of  Matrix- 
Isolated  HO^ 

The  theoretical  and  experimental  computer-simulated  spectra  of 

HC^  in  an  argon  matrix  [58]  were  obtained  according  to  the  procedure 

described  previously  for  the  computer  simulation  of  the  matrix 

spectrum  of  HOF.  The  computer-simulated  predicted  spectrum  of  HO^ 

is  shown  in  Fig.  3-6 (a).  Here  the  intensities  were  calculated  from 

the  ab  initio  4-31G  APTs  given  in  Table  3-21,  just  as  described 

previously  for  HOF.  Using  Eq.  (69)  we  calculated  the  (C&)  product 

needed  to  scale  the  plot  of  the  predicted  spectrum  by  fitting  the 

predicted  intensity  of  the  OH  stretching  mode  (3410  cm-1)  to  the 

observed  values  of  a and  Av  for  that  band  estimated  from 

mdx  1/2 

the  experimental  spectrum  of  HO^  in  the  argon  matrix  [58]  . Figure 
3 -6(b)  shows  the  computer  simulation  of  the  experimental  spectrum 
of  H02  in  the  argon  matrix.  The  predicted  band  intensity  for  the 
bending  anode  at  1390  cm  ^ is  within  factor  of  two  agreement  with 
the  corresponding  band  of  the  experimental  spectrum  and  the  under- 
estimated band  intensity  for  the  00  stretching  mode  is  shown  at 
1090  cm  . The  general  pattern  for  all  the  bands  is  consistent 
with  the  observed  spectrum  of  H02>  The  parameters  used  for  the 
computer  simulation  of  both  spectra  are  listed  in  Table  3-26. 

We  have  also  compared  the  theoretical  spectrum  predicted  from 
intensities  calculated  by  Komomicki  and  Jaffe  from  a TZPP  quality 
basis  set  [ 76  ] , with  the  experimental  matrix  spectrum  of  H02  in 
Fig.  3-7.  It  can  be  seen  that  even  this  extended  basis  set  does 
not  predict  very  well  the  general  pattern  for  all  the  bands  of  H02‘ 
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TABLE  3-26 

PARAMETERS  USED  FOR  THE  COMPUTER  SIMULATION  OF  THE  THEORETICAL  AND 
EXPERIMENTAL  SPECTRA  OF  MATRIX- ISOLATED  HO  (REFER  TO  FIG.  3-6) 


(C£=9.79 

Theoretical 
x 10-4  mmole  cm  ) 

(Cl= 

Experimental 
■9.79  x 10-4  irmole  cm. 

2> 

a 

V 

° -1 
(cm  X) 

AVl/2 

(orf-1) 

Q 

A (abs. ) 
(km  mole  "S 

v a 
(cm  ) 

v 

(cm~l) 

A° 

(km  mole  "S 

d 

a 

max 

3410.5 

6.0 

47.6 

3410.5 

6.0 

47.6 

0.73 

1389.0 

4.4 

78.5 

1389.0 

4.4 

148.3 

3.1 

1101.0 

4.0 

25.9 

1101.0 

4.0 

91.3 

2.1 

Wavenumber  position  of  maximum  absorption  from  the  experimental  matrix 

^spectrum  of  HO  as  reported  by  Smith  and  Andrews  in  ref.  [68] . 

Full  width  at  half-maximum  of  the  band  estimated  from  the  matrix- 

isolated  spectrum  [68] . 

^Taken  from  Table  3-25. 
d 

The  maximum  absorbances  of  the  HO  bands  estimated  from  the 
experimental  matrix  spectrum  of  HO^  in  the  argon  matrix  [68] . 
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In  particular,  the  relative  band  intensity  of  the  00  stretching 

is  also  underestimated.  The  parameters  used  for  the  computer 

simulation  of  both  spectra  are  shown  in  Table  3-27. 

Computer  Simulation  of  the  Matrix  Spectrum  of  H02  Predicted  by 
Transfer  of  APTs 

We  have  seen  previously  that  the  transferability  of  the  APTs 

depends  on  the  chemical  environment  of  the  atom  in  the  original 

molecule  and  that  in  the  molecule  to  which  the  APT  is  being  transferred. 

For  the  special  case  of  the  HO^  molecule,  the  best  approach  would  be 

to  transfer  the  APT  for  the  oxygen  atom  from  a system  containing  the 

00  bond,  and  the  APT  for  the  hydrogen  atom  from  H^O,  for  instance. 

Since  an  APT  for  a model  oxygen  atom  in  an  00  bond  has  not  yet  been 

determined  experimentally,  we  decided  to  investigate,  as  a first 

approximation,  whether  or  not  the  APTs  transferred  from  CH  OH,  (for 

3 

the  central  0 atom)  and  H^O  (for  the  H atom)  would  predict  successfully 
the  intensities  of  HO  . However,  to  predict  the  intensities  for  HO  by 
simply  transferring  the  APT  directly  from  HO  and  CH  OH  would  not 
be  realistic  since  the  central  oxygen  APT  is  expected  to  be  sub- 
stantially different  in  CH^OH  and  HO^ . To  prove  that  this  is  the 
case,  we  need  to  compare  the  contributions  to  the  ab  initio  APT 
for  the  oxygen  atom  in  both  molecules.  This  comparison  is  shown 
in  Table  3-28.  These  APTs  were  calculated  relative  to  the  coordinate 
system  shown  in  Fig.  3-8.  Notice  that  the  main  change  is  in  the 
charge  contribution  to  the  ab  initio  APT  of  the  oxygen  atom 
(-0.74  as  compared  to  -0.40  for  C^OH  and  HC>2,  respectively).  Notice 
also  that  the  charge-flux  and  overlap  APTs  are  different  for  the 
two  systems,  particularly  the  9p  (CF) , 9p  (CF)  and  9p  (0)  elements. 

X Z X 


9x 


9z 


9x 


Ch  O 
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TABLE  3-27 

PARAMETERS  USED  FOR  THE  COMPUTER  SIMULATION  OF  THE  THEORETICAL 
AND  EXPERIMENTAL  SPECTRA  OF  MATRIX- ISOLATED  HO  ( REFER  TO  FIG.  3-7  ) 


Theoretical  _2  Experimental 

(C$/=6.91  x 10“ 4 mmole  cm  ) (c£=6.91  x 10~4  mmole  cm  ) 


voa 
(cm  1) 

Avi/£  b 

(cm  ) 

Ac 

(km  mole 

V a 

, °-l. 
(cm  ) 

b 

(cm  1) 

AC 

(km  mole 

d 

a 

max 

) 

3410.5 

6.0 

67.5 

3410.5 

6.0 

67.5 

0.73 

1389.0 

4.4 

53.2 

1389.0 

4.4 

210.1 

3.1 

1101.0 

4.0 

46.  3 

1101.0 

4.0 

129.4 

2.1 

Wavenumber  position  of  maximum  absorption  from  the  experimental  matrix 
bspectrum  of  HO2  as  reported  by  Smith  and  Andrews  in  ref.  [68] . 

Full  width  at  half-maximum  of  the  band  estimated  from  the  matrix- 
isolated  spectrum  [68] . 

Taken  from  Table  3-25. 

The  maximum  absorbances  of  the  HO2  bands  estimated  from  the 
experimental  matrix  spectrum  of  HO2  in  the  argon  matrix  [68] . 
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TABLE  3-28 

CHARGE,  CHARGE-FLUX  AND  OVERLAP  CONTRIBUTIONS  TO  THE  AB  INITIO 
APT  ON  THE  CENTRAL  OXYGEN  ATOM  OF  CH  OH  AND  HO  (ATOM  1) 

(UNITS  ARE  e) 


— x 


CH3OHa 


p|0>  (C) 


(-0.19  0 

0 -0.07 

0.21  0 


(-0.86  0 

0 -0.76 

0.19  0 


(-0.40  0 

0 -0.40 

0 0 

(-0.19  0 

0 0 

0.02  0 

(0.003  0 

0 -0.04 

0.11  0 

(-0.59  0 

0 -0.45 

0.13  0 


-0.14  J 

0.12 

0 

0.39 

0.29 
0 

-0.15 


aTaken  from  Rogers  [35]  and  expressed  in  the  coordinate  system  shown 
^in  Fig.  3-8. 

cThis  work.  The  coordinate  system  is  shown  in  Fig.  3-8. 

Rogers  (private  communication)  and  expressed  in  the  coordinate 
system  shown  in  Fig.  3-8. 
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Coordinate  system  for  the  oxygen  atom  of  HOX,  X 


= CH  0. 

^ f 


Fig.  3-8. 
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However,  the  most  substantial  difference  is  indeed  in  the  charge 
contribution.  In  order  to  take  this  change  into  account  we 
scaled  the  experimental  APT  for  the  oxygen  atom  of  C^OH  according 
to 

P°(1)  (scaled)  = P°(1)  (CH  OH, Expl. )-  P0(1)  (CH  OH,  (C)  , ab  initio)  + 

^ X “X  -J 

P°  ^ (HO  , (C) ,ab  initio) 

X 2 

This  is  equivalent  to  adding  a correction  term  given  by 
P°^1)  (correction)  = p0^  (H  O,  (C)  ,ab  initio)  - 

X X 2 

(CH^OH,  (C)  ,ab  initio) 

to  the  experimental  APT  for  the  oxygen  atom  of  CH  OH.  We  have  used 

3 

the  same  procedure  to  transfer  the  hydroxy lie  hydrogen  APT  from  CH^OH 
and  the  oxygen  APT  scaled  also  from  CH^OH.  The  final  APTs  relative 
to  the  coordinate  system  shown  in  Fig.  3-4,  are  given  in  Table  3-29. 
The  predicted  intensities  for  these  transferred  APTs  are  shown  in 
Table  3—30  and  the  computer-simulated  spectra  in  Figs.  3-9  and  3-10. 
The  parameters  used  for  the  computer  simulation  of  both  spectra  are 
shown  in  Table  3-31. 

It  is  of  interest  now  to  compare  the  transferred  and  ab  initio 

APTs  for  H02,  which  is  also  presented  in  Table  3-29.  Notice  that 

there  are  differences  between  the  corresponding  elements  of  the 

hydrogen  APTs  and  oxygen  APTs  as  compared  to  the  corresponding 

elements  of  the  ab  initio  APTs.  The  difference  in  the  P°  ^ APT 
—x 

for  CH^OH  and  the  ab  initio  APT  is  due  to  the  fact  that  the  charge- 
flux  and  overlap  contributions  to  the  ab  initio  APT  of  CH^H  have 
not  been  corrected  (see  Table  3-28)  in  the  experimental  APT  for 
the  oxygen  atom  of  CH^OH.  Moreover,  there  are  differences  between 


TABLE  3-29 

TRANSFERRED  AND  AB  INITIO  4-31G  APTs  FOR  HO  (THE  COORDINATE  SYSTEM  IS  SHOWN  IN  FIG. 3-4  AND  UNITS 
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TABLE  3-30 

CALCULATED  INTENSITIES  FROM  TRANSFERRED  AND  AB  INITIO  APTS 
FOR  HO  (UNITS  ARE  km  mol-1) 


Band 

V (cm  1) 
k 

\a 

b 

\ 

c 

\ 

V (OH) 

3410.5 

20.2 

18.7 

47.6 

v2(A9) 

1389.0 

63.0 

30.3 

78.5 

v3 (00) 

1101.0 

12.6 

19.4 

25.9 

Calculated  intensities  by  transfer  of  the  APT  from  H^O  (for  the 
OH  hydrogen)  and  from  CH^OH  (for  the  OH  oxygen) . The  APTs  are 
presented  in  column  1 of  Table  3-29. 

Calculated  intensities  by  transfer  of  the  APT  from  CH^OH  (for  the 

OH  hydrogen  and  oxygen  atoms) . The  APTs  are  given  in  column  2 

of  Table  3-29. 
c 

Calculated  intensities  using  4-31G  APTs  given  in  column  3 of  Table 
3-29. 
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TABLE  3-31 

PARAMETERS  USED  FOR  THE  COMPUTER  SIMULATION  OF  THE 
MATRIX  SPECTRUM  OF  HO  PREDICTED  BY  TRANSFER  OF 

2APTs 


Predicted  by  transfer 

Predicted  by  transfer 

of  APTs 

for  H^O  and 

ch3oh 

of  APTs 

from  CH^OH 

(C&=9. 79 

x 10-4  mmole 

cm- 2) 

(C£=9.79 

x 10" 4 mmole 

cm  z) 

v a 

° -1, 
(cm  ) 

AVl/2 

(cm-1) 

c 

A 

(km  mole  "S 

v a 

, °-l. 
(cm  ) 

AVl/2b 

(cm-1) 

aC  -n. 

(km  mole  ) 

3410.5 

6.0 

47.6 

3410.5 

6.0 

47.6 

1389.0 

4.4 

63.0 

1389.0 

4.4 

30.3 

1101.0 

4.0 

12.6 

1101.0 

4.0 

19.4 

aWavenumber  position  of  maximum  absorption  from  the  experimental  matrix 
^spectrum  of  HO  as  reported  by  Smith  and  Andrews  in  ref.  [58] . 

°Full  width  at  half-maximum  of  the  band  [68]. 

Taken  from  Table  3-30. 
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( H) 

the  corresponding  elements  of  the  'APT  transferred  from  HO 

”” X 4 

or  CH^OH  as  compared  with  the  corresponding  elements  of  the  ab 
initio  APT  of  HO^ . Those  differences  reflect  in  the  intensity 
calculated  from  these  transferred  APTs.  For  instance,  Table  3-30 
shows  that  the  intensities  predicted  by  transfer  of  APT  from  H^O 
and  CH^OH  are  within  factor  of  two  or  better  agreement  with  the 
corresponding  predicted  intensity  from  ab_  initio  APTs,  except  for 
the  OH  stretching  mode . This  was  to  be  expected  due  to  the 
differences  between  the  APT  for  the  hydrogen  atom  of  H^O  and 
HO ^ (see  Table  3-29).  The  intensity  for  this  mode  is  within  factor 
of  2.4  agreement  with  the  corresponding  predicted  intensity  from 
ab  initio  APTs.  Table  3-30  also  shows  the  intensities  predicted  by 
transfer  of  APTs  from  CH^OH.  In  this  case,  the  intensities  predicted 
for  the  and  modes  of  HO^  are  within  factor  of  2.5  agreement  with 
the  corresponding  intensities  predicted  from  ab  initio  APTs . Only  the 
intensity  predicted  for  the  00  stretching  mode  is  in  good  agreement 
with  the  4-31G  intensity  (within  factor  of  1.3) . 

We  can  observe  that  the  general  pattern  of  the  spectrum  shown 
in  Fig.  3-9  is  very  similar  to  that  predicted  from  ab_  initio  APTs 
shown  in  Fig.  3-6.  Finally,  Fig.  3-10  shows  that  the  hydrogen  APT 
from  CH^OH  can  not  be  transferred  successfully  to  predict  the  in- 
tensities of  HO^.  We  shall  discuss  in  Chapter  4,  the  reasons  for 
nontransferability  of  the  hydrogen  APTs  in  terms  of  charge,  charge- 
flux  and  overlap  contributions  to  the  ab  initio  APT,  for  the  series 


of  molecules  studied  in  this  work. 
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Hydrogen  Hypochlorite 

As  was  previously  discussed  in  the  introduction  of  this 
work,  the  hydrogen  hypochlorite  molecule  is  one  of  the 
intermediate  species  in  the  cycle  of  the  chlorine  catalyzed 
decomposition  of  ozone  in  the  earth's  stratosphere.  We  have 
also  mentioned  that  the  controversy  about  the  photodissociation 
cross  section  of  HOC£  near  320  nm  has  recently  been  settled  by 
Molina,  Ishiwata  and  Molina  [13] . Their  studies  indicated 
that  HOC&  is  not  a long-lived  species  in  the  stratosphere 
and  therefore  it  does  not  play  an  important  role  as  a possible 
inert  reservoir  for  stratospheric  chlorine.  Nevertheless,  we 
think  that  our  ab  initio  SCF  calculation  of  the  intensity  of 
HOC£  will  be  of  interest  for  comparison  with  available  gas  phase 
experimental  intensities  [78]  and  with  the  ab  initio  SCF 
calculations  as  reported  by  Kcmornicki  and  Jaffe  [76].  Their 
calculation  [76]  was  carried  out  with  the  gradient  SCF  method 
using  a contracted  gaussian  basis  set  of  triple  zeta  accuracy  plus 
double  polarization.  We  have  performed  our  ab  initio  calculations 
using  a 4-31G  basis  set. 

Here  we  want  to  digress  for  a moment  to  point  out  that  we  have 
not  included  polarization  functions  in  the  basis  set  for  our 
calculations  in  HOC£  since  we  think  it  is  important  to  keep  the 
size  and  quality  of  the  basis  set  consistent  throughout  the  cal- 
culations if  meaningful  comparisons  are  to  be  made  of  the  predicted 
properties  for  a series  of  related  molecules.  Therefore,  we  want 


to  know  what  changes  in  the  contributions  to  the  ab  initio  APT  are 
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to  be  expected  when  the  chemical  environment  changes  within  a 
series  of  related  systems.  This  is  a point  which  we  want  to 
emphasize  throughout  this  study.  Only  by  calculations  of  the 
APTs  made  with  the  same  basis  set,  can  changes  in  these  properties 
be  investigated  and  analyzed.  For  instance,  it  is  not  appropriate 
to  calculate  the  intensity  for  one  particular  molecule  with  a 
minimum  ST0-3G  basis,  and  then  choose  an  extended  basis  for  the 
intensity  calculations  of  a second  molecule  and  a double-zeta 
quality  basis  set  for  a third  calculation  and  then  try  to  investigate 
the  variations  in  the  predicted  APTs  or  in  their  various  contributions 
when  the  chemical  environment  is  changed  fran  one  particular  molecule 
to  another.  We  also  want  to  point  out  that  the  inclusion  of 
polarization  functions  in  one  basis  set  may  improve  the  predicted 
intensity  for  one  specific  band;  for  example,  for  the  0-C#,  stretching 
mode.  However,  this  change  may  overestimate  or  underestimate  the 
intensities  for  other  bands  of  the  molecule,  depending  on  the 
particular  case  studied.  We  shall  discuss  this  matter,  which  is  of 
crucial  importance  in  more  detail  later. 

We  have  also  predicted  the  intensities  for  the  fundamental 
modes  of  HOC#,  by  transfer  of  the  APTs  from  other  molecules.  These  pre- 
dictions are  compared  with  the  ab  initio  and  experimental  results. 

These  studies  of  HOC#,  intensities  are  presented  in  the  next 
sections  of  this  chapter.  We  have  also  investigated  the  correlation 
of  the  ab  initio  charge-flux  contribution  to  the  APT  for  the  hydrogen 
atom  of  HOC#,  with  the  Mulliken  electronegativity,  as  shall  be 


discussed  in  Chapter  4. 
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Normal  Coordinate  Analysis 

The  geometrical  structure  of  HOC£  is  well  established  both 
by  microwave  [79,80]  and  infrared  spectroscopic  studies  [81] 
and  by  ab  initio  SCF  calculations  [11]  (see  Table  3-32).  In  our 
calculations  we  have  used  the  structural  parameters  given  by 

O 

Lindsey,  Lister  and  Millen  [79].  The  parameters  r(OC£)=1.69  A, 


r(OH)=0.97  A and  9 (H0c£)=103°  were  taken  from  reference[79] 
and  used  to  calculate  the  position  vector  for  the  HOC£  molecule 
as  given  in  Table  3-33,  according  to  the  coordinate  system  shown 
in  Fig.  3-11.  The  column  vector  given  by  Eq.  (2)  was  formed 
using  the  data  presented  in  Table  3-33.  This  column  vector, 
the  masses  of  the  atoms  (H=l. 007825,  0=15.994910  and  c£=34. 968854 
in  amu) , the  description  of  the  internal  coordinates  and  the  U 
matrix  were  used  as  input  to  program  WMAT  [27]  to  obtain  the 

and  G_  matrices,  In  Table  3-34  the  matrix  for  HOC£  is  presented. 
The  G matrix  was  then  diagonalized  by  program  WMAT  to  form  the 
W matrix,  according  to  Eq.  (28) . The  force  field  was  taken  from 
Schwager  and  Arkell  [82]  and  is  presented  in  Table  3-35.  This 
force  field  was  obtained  by  fitting  the  gas-phase  vibrational 
wavenumbers  of  HOC£  and  three  other  isotopically  substituted 
molecules,  using  the  least-squares  procedure  presented  in  Chapter  2. 
Table  3-35  also  compares  that  force  field  with  the  one  obtained 
by  Ogilvie  [50] . The  main  differences  between  those  force  fields 
is  observed  in  the  interaction  force  constants  F and  F „ 

r(OH),e  r'OCD.e 

We  decided  to  apply  both  force  fields  to  investigate  the  effect  in 


the  predicted  intensities.  We  find  no  appreciable  change  in  the 
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TABLE  3-32 

STRUCTURAL  PARAMETERS  OF  HOC£ 


r (OH)  a 

r (OC&) a 

0 (HOC&) 

References 

0.975 

1.689 

102.5° 

Mirri , Scappini  and 
Cazzoli , ref.  [80] 
(microwave) 

0.97 

1.69 

O 

103 

Lindsey,  Lister  and 
Millen,  ref.  [79] 
(microwave) 

0.97 

1.689 

104.8° 

Ashby  , ref.  [81] 
(infrared) 

0.97 

1.69 

O 

105 

Jaffe  and  Langhoff, 
ref.  [ 11  ] (ab  initi< 
SCF-CI  calculation) 

a . . 

Bond  length  m units  of  A 
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EQUILIBRIUM 

TABLE  3-33 

POSITION  VECTORS  FOR  HOC& (UNITS 

O 

ARE  A) 

Atom  No 

X 

Y 

z 

0,  1 

o 

• 

o 

o 

o 

o 

o 

H,  2 

o 

o 

-0.970000 

o 

#■ 

o 

Cl,  3 

0.0 

0.380167  -1 

.646685 

aThe  atom  numbering 

and  coordinate 

system  are  given 

in  Fig.  3-11. 
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Fig.  3-11.  Coordinate  system  and  atom  numbering  for  HOC£. 


TABLE  3-34 

THE  B MATRIX  FOR  HOC£ 
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TABLE  3-35 

THE  EXPERIMENTAL  F MATRICES  FOR  HOC& 


(UNITS  ARE  MDYNES  A-1) 


(1  aJ  A 2 = 1 mdyne  A“1  = 100  N m d) 


From  Schwagerc 
/ 7.: 


Vo. 


104 
0.0 
0 


From  Ogilvie 
/ 7.3177 
-0.1000 


\ 0.0 


0.0 

0.775 

0.677 


-0.1000 

0.7860 

0.3844 


0.0 

0.677 

3.980 


\ 

/ 


Taken  from  Schwager  and  Arkell  and  reported  in  ref.  [82] 

F11  = Fr(OH);  F22  = F0 ' F23  = Fr(OC£),0  F33  = Fr(OC£)  * 

^Taken  from  Ogilvie  and  reported  in  ref.  [50].  F, , = F 

11  r(OH) 

F12  = Fr(OH),0;  F22  = F0 ; F23  = Fr(OC2,),0  ^ F33  = Fr(OCil). 
Notice  that  the  following  force  constants  given  by  Ogilvie  are 
scaled  by  the  OCZ  and  OH  bond  lengths  of  HOC&,  according  to: 

F0  = F0/rlV  Frx0  = Fr;L0/ri  ^ Fr20  = Fr2e/r2'  where  rl  and  r2 


are  the  OH  and  OCi  bond  lengths  of  HOCit  respectively  and  the 
prime  indicates  here  the  scaled  force  constants  given  by  Ogilvie. 
Notice  also  that  the  force  constants  reported  by  Ogilvie  are  in 
units  of  N m--'-.  Here  F_  = F_,  since  the  U matrix  in  the  identity 
matrix. 
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predicted  intensities  for  HOC#.  that  can  be  attributed  to  the 
changes  in  the  force  field.  The  W matrix  obtained  from  program 
WMAT  and  the  F matrix  were  used  as  input  to  program  CHARLY  [27] 
to  obtain  the  eigenvectors  L_,  and  eigenvalues  A of  the  G F 
product  matrix.  The  13  matrix  generated  from  both  force  fields 
and  the  calculated  vibrational  wavenumbers  of  HOC#,  are  shown 
in  Tables  3-36  and  3-37,  respectively. 

Ab  Initio  4-31G  APTs  and  Calculated  Intensities  for  HOC#. 

The  APTs  for  HOC#,  were  quantum  mechanically  calculated  as 
before  for  HOF  and  HO^  using  a 4-31G  basis  set.  The  details  of 
this  calculation  have  been  discussed  in  Chapter  2.  The  APTs 
were  calculated  with  respect  to  the  coordinate  system  shown  in 
Fig.  3-11.  The  ab  initio  4-31G  APTs  for  HOC#,  are  given  in 
Table  3-38.  The  APT  for  the  central  O atom  (1)  was  calculated 
from  the  null  condition  given  by  Eq.  (62).  The  intensities  for 
each  fundamental  band  of  HOC#,  calculated  from  the  ab  initio 
4-31G  APTs  are  given  in  Table  3-39.  Also  shown  in  that  Table 
are  the  absolute  intensities  calculated  by  Komomicki  and  Jaffe  [76] 
with  the  GRADSCF  program  [37]  using  a contracted  gaussian  basis 
set  of  triple-zeta  accuracy  plus  double  polarization  (TZPP) . The 
experimental  intensities  for  HOC#,  gas  have  been  measured  by  Su 
et  al.  [78]  , with  a Fourier  transform  infrared  (FTIR)  spectrometer, 
and  these  values  are  also  included  in  Table  3-39.  This  Table 
shows  that  the  intensities  predicted  from  a 4-31G  basis  set  are 
within  a factor  of  two  agreement  with  the  corresponding  experimental 
intensities  except  for  the  OC#.  vibrational  mode.  The  inclusion  of 
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THE 

TABLE  3-36 

L MATRICES  FOR  HOC& (UNITS  ARE 

u'1/2 

1.026992 

0.006174 

0.001843 

■0.040833 

1.040859 

-0.291563 

-0.014374 

0.023306 

0.300611 

1.026785 

0.021614 

0.000500 

-0.057663 

1.074216 

-0.113130 

-0.013263 

-0.027549 

0.300303 

aL  matrix  generated  frcm  the  force  field  given  by  Schwager  and 
Arkell  [82] . 

L matrix  generated  frcm  the  force  field  given  by  Ogilvie  [50] . 


OBSERVED  AND  CALCULATED  VIBRATIONAL  WAVENUMBERS  FOR  HOC£ (UNITS  ARE  cm 
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TABLE  3-38 

AB  INITIO  APTs  FOR  HOC& (UNITS  ARE  e) 


— X 


p 


H 


— x 


PCA 

— x 


0 

-0.295 

0.209 

0 

0.244 

-0.105 

0 

0.0512 

-0.104 


Note:  A 4— 31G  basis  set  was  used  for  the  ab  initio  calculations 
of  the  APTs  for  HOC&.  The  P®  tensor  was  calculated  using 
the  null  condition.  The  coordinate  system  is  shown  in 
Fig.  3-11. 
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TABLE  3-39 


ABSOLUTE 

AND  EXPERIMENTAL 

INTENSITIES  OF 

HOC& (UNITS 

ARE  km  mole  ’*’) 

Band 

V (cm  'S 

A (4-31G) 3 

A (TZPP ) b 

A (Expl . ) C 

V (OH) 

3609.2 

75.3 

122.3 

56.6 

v (A0) 

1239.9 

39.3 

45.8 

73.9 

v (oc£) 

725.0 

0.089 

6.2 

10.6 

Calculated  intensities  using  the  4-31G  APTs  given  in  Table  3-38 
^(this  work) 

Calculated  by  Komomicki  and  Jaffe  and  reported  in  ref.  [ 76  ] using 

a contracted  gaussian  basis  set  of  triple— zeta  accuracy  plus  double 

polarization  (6s4p2d/5s3p2d/3s2p) , for  the  calculation  of  the  dipole 

moment  derivatives, 
c 

Experimental  absolute  intensities,  given  by  Su  et  al.  [78]. 
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polarization  function  is  definitely  important  for  predicting 
the  intensity  for  this  mode,  but  then  we  find  that  the  intensity 
for  the  OH  stretching  band  is  overestimated  as  has  been  previously 
discussed.  We  have  not  included  polarization  functions  in  our 
basis  set  for  reasons  also  discussed  in  the  introduction  to  the 
section. 

We  have  also  calculated  the  charge,  charge-flux  and  overlap 
contributions  to  the  4-31G  ab  initio  APTs  for  the  atoms  in  HOC&, 
according  to  the  CCFO  model  [52,53],  using  Eqs.  (63),  (64)  and  (67) 
respectively.  These  APTs  are  given  in  Tables  3-40,  3-41  and  3-42. 

The  coordinate  system  is  shown  in  Fig.  3-11.  In  Chapter  4,  we 
shall  compare  the  elements  of  those  APTs  for  all  of  the  molecules 
studied  in  this  work. 

Predicted  and  Experimental  Computer-Simulated  Gas-Phase  Spectra  of  HOCft 

The  theoretical  and  experimental  computer-simulated  spectra  of 
HOC£  in  gas  phase  are  shown  in  Figs.  3-12  and  3-13.  The  parameters 
used  for  the  computer  simulation  are  shown  in  Table  3-43.  Figure 
3-12  shows  the  predicted  spectrum  of  HOC£  from  ab  initio  4-31G  APTs 
compared  to  the  computer  fit  to  the  experimental  spectrum  of  HOC2, 
in  gas  phase.  Notice  that  the  intensity  predicted  for  the  OH 
stretching  at  3600  cm  ^ and  the  bending  mode  at  1240  cm  ^ are  within 
factor  of  two  or  better  agreement  with  the  corresponding  experimental 
intensities.  The  band  at  725  cm  ^ is  not  shown  in  that  spectrum 
since  its  intensity  is  only  0.089  km  mole  ^ as  predicted  from  the 


4-31G  APTs. 
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TABLE  3-40 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APTs 
ON  THE  HYDROGEN  ATOM  OF  HOC 5, 

(UNITS  ARE  e) 


PH(C) 
— x 


0.44 

0 

0 


0 

0.44 

0 


PH(CF) 
— x 


0 

0.12 

-0.06 


PH(0) 
— x 


0.04 

0 

0 


0 

-0.32 

-0.05 


ph(t) 

— x 


0.49 

0 

0 


0 

0.24 

-0.10 


Note:  The  charge  (C) , charge-flux  (CF)  and  overlap  (O)  contributions 
to  the  total  (T)  4-31G  APT  on  the  hydrogen  atom  were  calculated 
according  to  the  CCFO  model  [52,53].  The  coordinate  system  is 
shown  in  Fig.  3-11. 
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TABLE  3-41 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APTs  ON  THE 
CHLORINE  ATOM  OF  HOC& (UNITS  ARE  e) 


pC£ 
— x 


(C) 


0 

0.22 

0 


nO 

P (CF) 
— x 


0 

-0.02 

0.09 


pa 

— x 


(O) 


0 


-0.15 


-0.19 


P<*(T> 

— X 


0.07 


0 

0.05 


-0.10 


-0. 19 


Note:  The  charge  (C) , charge-flux  (CF)  and  overlap  (O)  contributions 
to  the  total  (T)  APT  on  the  chlorine  atcm  were  calculated 
according  to  the  CCFO  model  [52,53].  The  coordinate  system 
is  shown  in  Fig.  3-11. 
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TABLE  3-42 

CHARGE,  CHARGE-FLUX,  OVERLAP  AND  TOTAL  APTs  ON  THE 
OXYGEN  ATOM  OF  HOCil  (UNITS  ARE  e) 


P°  (C) 
— x 


0 

P (CF) 
—x 


0 

-0.10 

-0.03 


P°  (0) 
—x 


0 

0.47 

0.24 


P°  (T) 
-x 


0 

-0.30 

0.21 


Note:  The  charge  (C) , charge-flux  (CF)  and  overlap  (O)  contributions 
to  the  total  (T)  4-31G  APT  on  the  oxygen  atom  were  calculated 
according  to  the  CCFO  model  [52,53].  The  coordinate  system  is 
shown  in  Fig.  3-11. 


Fig.  3-12.  Computer  simulation  of  the  predicted  spectrum 
of  HOC 5-  from  4-31G  APTs  and  the  computer- 
simulated  gas  phase  spectrum  of  HOC&. 

(a)  Predicted 

(b)  Experimental 
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Fig.  3-13.  Computer  simulation  of  the  predicted  spectrum 
of  HOCii  using  a contracted  gaussian  basis  set 
of  triple-zeta  accuracy  plus  double  polariza- 
tion and  the  computer-simulated  gas  phase 
infrared  spectrum  of  HOC£. 

(a)  Predicted  (Komomicki  and  Jaffe  in  ref. 

[76]). 

(b)  Experimental. 
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TABLE  3-43 

PARAMETERS  USED  FOR  THE  COMPUTER  SIMULATION  OF  THE  PREDICTED  AND  EXPERIMENTAL 

GAS  PHASE  SPECTRA  OF  HOC& 
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Figure  3-13  shows  the  computer  simulation  of  the  infrared 
spectrum  of  HOCil  calculated  from  the  intensities  reported  by 
Komornicki  and  Jaffe  [76] . Notice  that  the  OCi  stretching  mode 
now  appears  with  appreciable  predicted  intensity  at  725  cm  1. 

The  bands  at  1240  and  725  cm  ^ are  in  very  good  agreement  with 
the  experimental  intensities.  However  the  intensity  for  the  OH 
stretching  mode  is  overestimated.  Kanornicki  and  Jaffe  suggested 
that  electron  correlation  effects  could  make  a large  contribution 
to  the  calculated  dipole  derivatives  for  the  OH  stretching 
mode  in  HOCil. 

It  is  well  known  that  calculations  at  the  SCF-CI  level  with 
an  extended  basis  set  including  polarization  functions  result  in 
reduced  Hartree-Fock  force  constants  and  dipole  moments  [83] . 
However,  it  is  also  true  that  calculations  of  dipole  moments 
derivatives  are  very  sensitive  to  the  quality  of  the  basis  set  used. 
Let  us  consider,  for  example,  the  predicted  intensity  for  the 
vibrational  mode  of  H^O.  This  is  shown  in  Table  3-44,  where  the 
predicted  intensity  for  this  mode  are  ordered  according  to  the 
quality  of  the  basis  set.  Notice  that  1.07  km  mole  ^ is  the 
intensity  predicted  with  a 4-31G  basis  set.  The  next  three  sub- 
sequent intensity  values  were  predicted  with  very  large  basis  sets 
of  increasing  quality  at  Hartree-Fock  level  and  including  different 
sets  of  polarization  functions  [84-86]  . 

That  Table  also  includes  several  experimental  values  reported 
for  the  intensity  of  this  mode  [87-89] . It  is  interesting  now 
to  notice  that  by  increasing  the  quality  of  the  basis  set  and 
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TABLE  3-44 

COMPARISON  OF  THE  SIZE  AND  QUALITY  OF  A BASIS  SET  WITH  THE 


PREDICTED  INTENSITY 

FOR  THE  V 

VIBRATIONAL  MODE  OF  ^0 

VH20)* 

E 

o 

Description 

(km  mole  ) 

(Hartrees) 

1.07a 

-75.909 

SCF  level (4- 3 1G  basis) 

8.77b 

-76.035 

SCF-HF  level 

13. 6C 

-76.041 

SCF-HF  level 

15.  4d 

-76.051 

SCF-HF  level 

(2.24)®(2.51)5  (2.26)9 

-76.431 

Experimental 

aQuoted  by  Zilles  [5:3,p.73],  The  intensities  for  the  and 
modes  are  within  factor  of  two  agreement  with  the  corresponding 
experimental  values. 

Pul ay  [84] . 

d Smith,  Jorgensen  and  Ohm  [85], 

Krohn  and  Kern  [ 86] . 

^Clough  et  al.  [87], 

Toth  [ 88] . 

^Flaud  and  Camy-Peyret  [89]. 
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performing  the  calculations  at  the  Hartree-Fock  level  causes  the 
intensities  to  be  overestimated  by  factors  of  4,6  and  7.  However, 
the  predicted  intensity  with  the  4-31G  basis  is  within  factor 
of  two  agreement  with  the  experimental  results. 

This  does  not  seem  to  be  attributable  to  the  error  in  the 
Hartree-Fock  approximation  (dissociation  to  ions  instead  of  neutral 
atoms) . Instead,  we  believe  this  error  is  due  to  the  criterion 
chosen  to  measure  the  quality  of  the  basis  set:  namely,  that  it 
produces  a lower  energy  eigenvalue.  It  is  unrealistic  to  believe 
that  basis  sets  including  different  sets  of  polarization  functions 
which  have  been  optimized  only  from  energy  eigenvalue  considerations 
would  also  necessarily  reflect  the  proper  flexibility  to  accommodate 
the  electronic  charge-flux  redistribution  when  an  atom  is  displaced 
from  its  equilibrium  configuration. 

In  fact,  we  suggest  that  the  discrepancies  in  the  predicted 
values  for  the  dipole  moment  derivatives  with  different  basis  sets 
(for  example  those  from  a double  or  triple-zeta  plus  polarization 
basis  set  compared  with  results  from  a 4-31G  basis  for  H^O,  HOC& 
and  HC^ ) can  be  understood  on  these  terms.  The  inclusion  of 
polarization  functions  can  be  very  important  for  lowering  the  total 
molecular  energy  but  it  does  not  assure  the  use  of  a perfect 
formally  balanced  basis  set  for  dipole  moment  derivative  calculations. 
In  other  words , the  orbital  exponents  of  the  polarization  functions 
are  not  scaled  to  reflect  the  proper  charge-flux  redistribution  in 
a molecule  for  a bond  stretching  or  bending  deformation.  It  is 
important  to  point  out  that  energy  eigenvalues  (variation  theorem) 
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should  not  be  used,  in  general,  as  the  sole  criterion  for  quantum 
mechanical  calculations  of  properties  depending  as  sensitivily  as 
dipole  moments  on  the  charge  distribution. 

Mulliken  [55]  was  the  first  to  recognize  the  difficulty  of 
choosing  a formally  balanced  polarization  set  for  heteiopolar 
molecules  and  has  pointed  out  that  there  is  no  assurance  that 
2p  polarization  function  on  the  hydrogen  of  HF  would  provide  a 
perfect  balance  for  the  d polarization  set  of  the  fluorine  atom. 
Clementi  [90]  has  also  pointed  out  that  if  the  dipole  moment  would 
have  been  included  as  an  additional  criterion  for  the  optimization 
of  the  orbital  exponents  of  the  basis  set  chosen,  a better  agreement 
between  experiment  and  theory  could  have  been  obtained  for  the 
dipole  moment  calculation  of  HF.  However,  it  is  clear  that  dipole 
moment  can  not  be  the  best  criterion  for  choosing  a quality  basis 
set,  since  it  is  not  an  atomic  property  and  therefore  different 
sets  of  orbital  exponents  would  have  to  be  used  for  different 
molecules.  However,  it  is  possible  that  an  extended  flexibility 
of  a basis  set  could  be  achieved  at  the  SCF-HF  level  if  intensity 
parameters  could  be  used  as  an  additional  criterion  in  choosing 
a quality  basis  set.  Accordingly,  we  suggest  that  effective 
atomic  charges  could  be  used  as  an  additional  criterion  for  optimiza- 
tion of  the  outer-shell  orbital  exponents  of  a basis  set,  since 
this  atomic  property  is  invariant  to  the  normal  coordinate  analysis 
and  seems  to  be  transferable  among  a great  number  of  molecules  [52] . 

In  Chapter  4 we  investigate  further  the  prediction  of  the 
spectrum  of  HOCfi.  using  a 4-31G  basis  set,  without  the  inclusion  of 
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polarization  functions.  It  is  shown  that  if  a reasonable 
correction  is  made  to  the  charge-flux  tensor  of  the  hydrogen 
atom,  the  predicted  spectrum  almost  superimposes  on  the 
corresponding  experimental  spectrum  of  HOC£. 

Computer  Simulation  of  the  Gas-Phase  Spectrum  of  HOCJl  Predicted 
by  Transfer  of  APTs 

In  this  section  we  shall  discuss  several  attempts  to  transfer 
experimental  APTs  from  related  molecules  to  predict  the  spetrum 
of  HOCJl.  This  study  involves  APTs  transferred  from  four  sets  of 
molecules:  (C  Z CN  [91],  HO  [53]),(CH  C£  [34],  HO  [53]  ) , (C  % CN  [91]  , 
CH3OH  [35]),  and  (CH3C  i [34],  CH  OH  [35]).  The  APTs  transferred 
from  those  molecules  are  shown  in  Tables  3-45,  3-46,  3,47  and  3-48 
respectively,  and  the  coordinate  system  relative  to  those  APTs  is 
shown  in  Fig.  3-11.  The  predicted  intensities  are  compared  in 
Table  3-49  and  the  computer  simulations  of  the  predicted  and  ex- 
perimental spectra  are  shown  in  Figs.  3-14  to  3-17.  The  parameters 
used  for  the  computer  simulation  are  listed  in  Table  3-50. 

It  is  interesting  to  compare  the  elements  of  the  chlorine  APTs 
for  C£CN  and  CH^Cil.  Notice,  for  instance,  that  the  elements  of  the 
chlorine  APTs  for  C&CN  and  CH^CJl  as  given,  for  instance  in  Tables 
3-45  and  3-46  are  substantially  different.  These  differences  are 
expected,  since  a chlorine  atom  bonded  to  the  CN  group  is  chemically 
quite  different  from  one  bonded  to  the  methyl  group.  The  most 
substantial  differences  are  observed  for  the  3p^/3x  and  3p^/3y 
elements.  Notice  also,  that  the  signs  of  the  derivatives  for  those 
elements  are  also  different.  It  is  also  interesting  to  notice  that 
the  Sp^/Sz  elements  are  approximately  the  same  for  both  molecules. 
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TABLE  3-45 

TRANSFERRED  EXPERIMENTAL  APTs  FROM  CilCN  AND  H 0 

(UNITS  ARE  e)  2 


P° 

x 


PH 

—x 


pC£ 

—x 


0 

-0.246 

-0.078 

0 

0.132 

-0.049 

0 

0.114 

0.127 


Note:  P^  was  calculated  from  the  null  condition  given  by  Eq.  (62) . 

PH  was  transferred  from  H„0  [53] . 

— x 2 

pC&was  transferred  from  C&CN  [ 91  ] . The  original  APTs  were 
transformed  to  the  coordinate  system  shown  in  Fig.  3-11. 
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TABLE  3-46 

TRANSFERRED  EXPERIMENTAL  APTs  FROM  CH  C£  AND  H 0 

(UNITS  ARE  e)  3 2 


P° 

-x 


■0.096 

0 

0 


PH 

-x 


c£ 

P 

—x 


0 

0.113 

0.002 

0 

0.132 

-0.049 

0 

-0.245 

0.047 


Note:  P^  was  calculated  from  the  null  condition  given  by  Eq.  (62) . 

PH  was  transferred  from  HO  [53] . 

~~ 2 

P was  transferred  frcm  CH  C£[34]. 

T^ie  coordinate  system  is  shown  in  Fig.  3-11. 
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TABLE  3-47 

TRANSFERRED  EXPERIMENTAL  APTs  FROM  C^CN  AND  CH  OH 

(UNITS  ARE  e) 


P° 

—x 


PH 

—x 


0 

-0.255 

-0.145 

0 

0.141 

0.018 


PCi 

—x 


0 

0.114 

0.127 


Note:  P was  calculated  from  the  null  condition  given  by  Eq.  (62). 
— ^ 

P was  transferred  from  CH  OH  [35] . 

3 

P was  transferred  frcm  C&CN  [91] . 

T^ie  coordinate  system  is  shown  in  Fig.  3-11. 
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TABLE  3-48 

TRANSFERRED  EXPERIMENTAL  APTs  FROM  CH  C£  AND  CH  OH 
(UNITS  ARE  e) 


— x 


c£ 

p 

— x 


-0.026 

0 

0 

0.260 

0 

0 

-0.234 

0 

0 


0 

0.104 

-0.065 

0 

0.141 

0.018 

0 

-0.245 

0.047 


Note:  P was  calculated  from  the  null  condition  given  by  Eq. 
— ^ 

P was  transferred  from  CH_0H  [35] . 

3 

P was  transferred  from  CH^C#-  [34]  . 

T?ie  coordinated  system  is  shown  in  Fig.  3-11. 


(62)  . 


TABLE  3-49 

PREDICTED  INTENSITIES  FROM  TRANSFERRED  APTs  FOR  HOC£ (UNITS  ARE  km  mole 
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Fig.  3-14. 


Computer  simulation  of  the  predicted  spectrum 
of  HOC£  from  transferred  APTs  from  C£CN  and  H^O 
and  the  computer- simulated  gas  phase  spectrum 
of  HOC£. 

(a)  Predicted 


(b)  Experimental 
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Fig.  3-15. 


Computer  simulation  of  the  predicted  spectrum 
of  HOCZ  from  transferred  APTs  from  CH^C&  and 
HO  and  the  computer- simulated  gas  phase  spectrum 
of  HOC£. 

(a)  Predicted 


(b)  Experimental 
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Fig.  3-16.  Computer  simulation  of  the  predicted  spectrum 
of  HOC 5.  from  transferred  APTs  from  CilCN  and 
CH^OH  and  the  computer- simulated  gas  phase 
spectrum  of  HOC#.. 

(a)  Predicted 

(b)  Experimental 
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Fig.  3-17.  Computer  simulation  of  the  predicted  spectrum 
of  HOCS/  from  transferred  APTs  from  CH^CJl  and 
CH^OH  and  the  computer- simulated  gas  phase 
spectrum  of  HOCil. 

(a)  Predicted 

(b)  Experimental 


ANSMI TTANCE  TRANSMITTANCE 

•25  0.50  0.75  1.00  JJ.00  0.25  0.50  0.75  1.00 
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FREQUENCY  (CM  (-11 ) (X10l  ) 


CC  o' 


H 1- 


-t- 


-f- 


— (- 

ms.  oo 


— f- 

102.50 


-+- 


^100.00  357.50  315.00  272.50  230.00  187.50 

FREQUENCY  (CM  (-1) ) (X 1 0l  ) 


6( 


TABLE  3-50 

PARAMETERS  USED  FOR  THE  COMPUTER  SIMULATION  OF  THE  GAS-PHASE  SPECTRUM  OF  HOC£ 

PREDICTED  BY  TRANSFER  OF  APTs 


142 


E 

0 

i — i 

1 

oC 


(0 
P> 

c 

03 
E 

•H 
P 
03 

ft  I 

X 

H 


3 

M 


C 


CN  i— I 

in  'e 

P>  u 
< — 


(0  i 

;>  u 


o 


< 


33 

O 

sT 

u 

O 

ft 

33 

U 


33 

O 

c 

as 

u 

z 

u 

u 


CN 


0) 

P> 

o 


03  *1 
03  I 
P O 
ft  i— I 


U 


as 

u 

c* 

X 

o 


CN 


as 

z 

u 

o 


CN  i 


< 


E 

o o 


VO 

<T» 

VO 

• 

c 

VO 

00 

o 

• 

•H 

LD 

rH 

r— i 

o 

'0 

in 

a; 

O 

o 

o 

■ — > 

pi 

• 

• 

• 

p 

O 

o 

o 

CJ 

o 

r- 

LO 

LO 

o 

ft 

SB 

a) 

p 

4H 

o 

u 

CN 

<T> 

o 

s 

o 

• 

• 

• 

p 

as 

CT» 

LO 

P 

o 

00 

CN 

P> 

4a 

VO 

CN 

O 

0 

ro 

rH 

Q) 

ft 

£ 

03 

g 

p 

03 

pi 

03 

o 

A 

03 

aS 

ft 

CO 

00 

00 

ft 

03 

• 

• 

• 

£0 

00 

r- 

03 

03 

rH 

00 

03 

03 

O' 

03 

as 

i— 1 

ft 

(0 

pi 

03 

c 

03 

03 

ft 

E 

*H 

03 

P 

as 

0) 

pi 

ft 

r— 1 

vO 

co 

X 

E 

• 

• 

• 

<13 

0 

00 

rH 

CD 

p 

CN 

00 

<D 

4a 

as 

P> 

03 

03 

E 

P> 

0 

03 

p 

E 

• 

4a 

*H 

pj 

ft 

c 

03 

0 

03 

i 

CN 

CN 

ft 

•rH 

n 

• 

• 

• 

+J 

03 

VO 

V0 

00 

ft 

C 

03 

rH 

^0 

p 

0) 

i— 1 

0 

aa 

A 

03 

<3 

aa 

03 

Ea 

03 

as 

P3 

E 

g 

0 

3 

4a 

P 

VO 

in 

£ 

0 

4a 

• 

• 

• 

•rH 

rH 

LO 

00 

X 

£ 

C 

CN 

VO 

03 

3 

03 

E 

E 

A 

•H 

<3 

pi 

X 

pi 

0 

03 

— 

E 

C 

i 

i — i 

0 

4a 

i 

•H 

r— l 

03 

O 

o 

o 

pi 

o5 

i — 1 

• 

• 

•H 

as 

Q 

O 

o 

o 

cn 

E 

r- 

m 

in 

0 

pi 

fa 

o3 

E 

.X 

u 

as 

• 

CD 

p 

i— . 03 

Q 

03 

00  P 

•H 

r~  oj 

p 

5 

1 — * 

CN 

<S\ 

o 

c 

03 

• 

• 

• 

Q) 

rH 

• P> 

(7> 

(T> 

m 

> 

rH 

4a  -H 

O 

oo 

CN 

3 

03  a 

vO 

CN 

£ 

fa 

P D 

00 

rH 

<o  aa 

o 

143 


The  absolute  difference  is  only  0.03e,  which  is  within  the  exper- 
imental error. 

From  our  previous  discussion  about  the  transferability  of  APTs , 
we  expect  that  the  predicted  intensities  for  HOCil  from  APTs  trans- 
ferred from  the  (CilCN,  H^O)  pair  may  be  in  better  agreement  with 
experiment  than  predictions  from  any  of  the  other  sets  of  trans- 
ferred APTs.  This  comparison  is  shown  in  Table  3-49.  Notice  that 
the  predicted  intensities  are  in  general  within  factor  of  two 
agreement  with  the  experimental  intensities  for  intensities  pre- 
dicted from  CilCN  and  H^O  (see  also  Fig.  3-14).  However,  in  spite 
of  our  arguments  concerning  expected  transferability,  the  results 
in  Table  3-49  for  all  four  sets  of  parameters  are  close  enough 
to  the  experimental  values  that  we  can  not  choose  between  them. 


CHAPTER  4 

AN  ANALYSIS  OF  THE  ATOMIC  POLAR  TENSORS 
Introduction 

In  the  last  chapter  we  have  examined  the  differences  between 
the  elements  of  the  APTs  for  several  molecules , to  find  that 
some  APTs  are  transferable  and  others  are  not,  depending  upon  the 
chemical  environment  of  the  atom  or  group  of  atoms  from  which  the 
APT  is  being  transferred  and  that  of  the  new  molecule.  We  have 
been  able,  therefore,  to  predict  successfully  the  relative  or 
absolute  intensities  of  molecules  such  as  H00  and  HOC&,  by  proper 
transfer  of  the  APTs  from  ^0,  CH^OH  and  C&CN.  This  last  example  shows 
that  APTs  are  indeed  transferable  intensity  parameters,  providing 
that  some  judgement  is  used  in  the  transfer. 

We  have  also  shown  that  our  ab  initio  calculations  using  a 
4-31G  basis  set  gave  predicted  absolute  intensities  that  are 
accurate  within  the  experimental  uncertainties , as  do  also  the 
predictions  made  by  transferring  APTs.  We  have  found  that  only 
the  intensity  of  the  OCil  stretching  mode  for  H0C£  is  poorly  predicted 
with  the  4-3 1G  basis.  We  have  also  compared  the  predicted  spectrum 
for  H0C£  using  the  4-31G  basis  set  with  that  predicted  [76]  with 
a triple-zeta  accuracy  basis  set  plus  double  polarization.  In  this 
chapter  we  are  interested  in  investigating  the  changes  in  the 
fluorine  APT  of  the  CH^F  and  HOF  molecules  which  explain  why  this 
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tensor  is  not  transferable,  as  well  as  the  changes  in  the  hydrogen 

APTs  for  the  HOX  series  of  molecules  which  take  place  as  a result 

of  the  perturbing  X atom  in  terms  of  the  CCFO  quantum  mechanical 

model  for  the  APT  discussed  in  the  last  chapter.  As  a result  of 

the  HOX  study,  we  shall  find  a correlation  between  the  overlap 
H H 

(P  (0) ) and  charge-flux  (P  (CF))  contributions  to  the  APT  for  the 
—x  —x 

hydrogen  atom  and  the  Mulliken  electronegativity  of  the  X atom. 

Using  this  correlation  to  modify  the  APT  calculated  for  the 

hydrogen  atcm  in  HOC£,it  was  possible  to  predict  accurately  the 

spectrum  of  HOC£.  The  overlap  tensor  from  the  4-31G  calculation 

of  the  hydrogen  atom  in  HOC&  already  obeys  this  correlation.  The 

charge  tensor-contribution,  P (C) , to  the  APT  for  the  hydrogen 

—x 

atom  was  found  to  be  almost  constant  regardless  of  the  X atom, 
except  where  X=H  for  the  H^O  molecule. 

We  believe  that  the  4-31G  basis  does  not  accurately  predict 
the  polar  tensor  for  the  atoms  in  HOCil  because  the  basis  set  should 
include  polarization  functions.  Hence  it  seems  reasonable  to  use 
the  correlations  with  electronegativity  established  for  the  more 
accurately  calculated  APTs  for  HOX  molecules  to  correct  the  APT 
for  the  hydrogen  atom.  We  think  that  it  would  be  possible  to 
improve  the  basis  set  and  make  more  accurate  quantum  mechanical 
calculations  by  a simple  modification  of  the  basis  set  used. 

Namely,  we  should  optimize  the  orbital  exponents  of  added  polariza- 
tion functions  with  respect  to  some  intensity  parameter  (here 
we  suggest  the  atomic  effective  charges)  and  then  use  a standard 
molecular  set  of  these  scaled  factors  based  on  an  average  of  the 
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optimum  values  obtained  from  a series  of  model  molecules. 

This  additional  criterion  could  prove  to  be  extremely  important 
for  calculations  of  dipole  moments  and  derivatives,  within  the 
Hartree-Fock  formalism.  This  new  basis  set  will  reflect  criteria 
based  both  on  energy  eigenvalue  calculations  and  also  on  dipole 
moment  derivative  calculations. 

In  the  next  section  we  shall  compare  the  contributions  of 
the  CCFO  terms  to  the  ab  initio  4-31G  APTs  for  CH^F  and  HOF. 
Following  that  we  show  how  we  can  proceed  to  modify  the  charge- 
flux  contributions  to  the  hydrogen  atom  APT  for  the  HOX  series 
of  molecules,  based  on  these  electronegativity  correlations. 

Comparison  of  the  Ab  Initio  4-3 1G  APTs  for  the 

Fluorine  Atom  in CH-^F  with  that  in  the  HOF  Molecule 

We  have  seen  in  the  previous  chapter  that  the  fluorine  APT 
from  CH^F  can  not  be  successfully  transferred  to  predict  the 
intensities  for  HOF  (see, for  example,  Table  3-14).  We  have  shown 
that  the  predicted  intensities  for  HOF  by  transfer  of  the  fluorine 
APT  from  CH^F  and  the  hydrogen  APT  from  H^O  are  in  poor  agreement 
with  intensities  predicted  from  ab  initio  APTs  for  HOF  or  from 
relative  experimental  intensities  of  this  molecule.  Accordingly, 
the  intensities  predicted  from  APTs  tranferred  frcm  CH^F  and  H^O 
show  an  order  of  increasing  intensities  frcm  high  wavenumber 
(OH)  to  low  (OH)  precisely  the  opposite  of  what  is  observed 
experimentally  or  predicted  quantum  mechanically.  The  reason  for 
the  nontransferability  of  the  fluorine  APT  from  CH^F  to  HOF 
shows  in  the  contributions  to  the  ab  initio  APT  for  both  molecules , 
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as  compared  in  Table  4-1.  This  Table  shows  that  the  charge- 
tensor  contributions  are  substantially  different  for  the  two 
F atoms.  The  fluorine  charge  tensor  gives  the  contribution 
from  the  gross  atomic  charge  on  the  fluorine  atcm  at  equilibrium 
molecular  configuration  and  reflects  the  chemical  differences 
between  the  F atoms  jn. the  CH^F  and  HOF  molecules.  Therefore, 
the  fluorine  APT  of  CH^F  and  HOF  is  different  in  each  molecule 
and  obviously  not  transferable.  Notice  that  the  diagonal 
elements  9p  of  the  charge-flux  and  overlap  tensor  contributions 


itself  an  answer  for  those  who  still  believe  that  the  fluorine 
atom  is  so  electronegative  that  it  does  not  matter  in  its  electronic 
properites  just  what  atom  it  is  attached  to. 


Before  considering  in  some  detail  the  electronegativity 
correlation  with  the  charge-flux  tensor  of  the  hydrogen  atom  for 
the  HOX  series  of  molecules  we  wish  to  point  out  the  logical  reason- 
ning  which  led  us  to  this  study. 

We  have  previously  noticed  that  the  experimental  chlorine  APT 
for  C5.CN  and  the  ab  initio  predicted  APT  for  the  chlorine  atom  in 
H0C5,  are  quite  different.  Therefore  it  would  be  of  interest,  for 
example,  to  compare  both  APTs  with  the  ab  initio  APT  of  the  chlorine 
atom  in  the  diatomic  0C5.  molecule.  This  comparison  is  shown  in 
Table  4-2 . The  APT  for  the  chlorine  atom  in  diatonic  0C5, 


Predicted  Vibrational  Intensity  and  the  APTs 
For  Diatomic  0C5, 
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TABLE  4-1 

COMPARISON  OF  THE  CHARGE,  CHARGE -FLUX  AND  OVERLAP 
CONTRIBUTIONS  TO  THE  AB  INITIO  4-31G  APT  FOR  CH  F AND  HOF 


— x 


CH3F 


HOFC 


(-0.46  0 

0 -0.46 

0 0 

(0.02  0 

0 -0.24 

0 0 


PF(0) 

-x 


0 0 
-0.35  -0.04 


0.01 


0.02 


F 

P (T) 
— x 


/- 


0.38 

0 


0 

-0.87 


0 

0 


0 


-0.38 


0 

-0.44 

0.05 


0 \ 

-0.16 

-0.17  / 


aThe  charge  (C) , charge-flux  (CF)  and  overlap  (O) , contributions 
to  the  total  4-31G  APT  on  the  fluorine  atom  were  calculated 
according  to  the  CCFO  model  [52,53].  The  coordinate  system  is 
^shown  in  Fig.  3-1. 

The  fluorine  APT  from  CH^F  was  from  Chin  (private  communication 
of  work  in  progress  at  the  University  of  Florida) . 

This  work. 
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TABLE  4-2 

COMPARISON  OF  THE  CHLORINE  APT  FOR  SEVERAL  MOLECULES 
(THE  COORDINATE  SYSTEM  IS  SHOWN  IN  FIG. 3-11) 


p<CA) 

— x 


(0.164  0 0 

0 0.133  0.135 

0 0.135  -0.423 


(0.143  0 0 

0 0.114  0.127 

0 0.127  -0.407 


(0.0743  0 0 

0 0.0512  0.021 

0 -0.104  -0.190 


(-0.234  0 0 

0 -0.245  0.047 

0 0.047  -0.435 


oc£a 


c£cn 


b 


hoc£.c 


CH3c£d 


aAb  initio  4-31G  APT  (this  work). 

Experimental  chlorine  APT  reported  in  ref.  [ 91]  and  expressed 
here  in  the  coordinate  system  shown  in  Fig.  3-11 (see  also 
Table  3-45  of  this  work) . 

Ab  initio  4-31G  APT  taken  from  Table  3-38  (this  work) 
dTaken  from  Table  3-46. 
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can  be  written  [1,92] . 

') 

3pz/8z  I 

Here  p is  the  magnitude  of  the  permanent  dipole  moment  vector 

at  equilibrium  molecular  configuration,  R is  the  equilibrium 

e 

intemuclear  distance  in  0C£  and  3p  /3  is  the  magnitude  of  the 

z z 

dipole  derivative  with  respect  to  the  z axis.  The  experimental 

values  for  p°  and  R were  obtained  from  ref erence [93 ] 

e 

[p°=1.239  D (compared  to  a calculated  value  of  p°=1.112  D)  and 

O 

R =1.569  A]  and  the  3p  /3z  element  of  this  tensor  was  calculated 
e z 

as  described  for  the  molecules  using  a 4-31G  basis  set.  Table  4-2 
shows  the  comparison  of  this  APT  for  Ci  in  0C£,  C&CN,H0C&  and  CH^Cii. 
The  coordinate  system  is  shown  in  Fig.  3-11.  Notice,  for  instance, 
that  the  chlorine  APTs  from  OCJl  and  C&CN  are  quite  similar,  in 
contrast  with  the  predicted  chlorine  APT  for  H0C&.  It  seems  clear 
that  the  chlorine  APT  of  H0C£  is  not  well-predicted  using  a 4-31G 
basis  set,  probably  due  to  the  lack  of  polarization  functions  on  the 
hydrogen  and  chlorine  atoms. 

Here  it  is  interesting  to  mention  that  the  Pauling  electro- 
negativity of  the  0 atom  and  the  electronegativity  of  the  CN  group 
are  practically  the  same  (3.44  compared  to  3.3  [94])  and  the  polarity 
of  the  C C&  bond  in  both  molecules  is  the  same  (c£+0  , C£+CN  ) . 

For  the  CH^CJl  molecule  the  polarity  of  the  C C£  bond  (C+C£  ) is 
the  opposite  to  that  of  C£o  or  CitCN.  However,  we  have  not  observed 
in  Figs.  3-14  and  3-15  any  substantial  effect  on  the  predicted 
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spectrum  of  HOC&  from  APTs  transferred  from  either  molecule. 
Nevertheless,  we  think  the  C&  atom  in  HOC&  is  expected  to  be 
more  like  those  in  C£CN  or  cHo  than  like  that  in  CH^CJi. 

In  the  next  section  we  shall  discuss  how  we  modified  the 
charge-flux  tensor  of  the  hydrogen  atom  in  HOC&  to  predict  the 
intensities  for  this  molecule,  using  also  the  chlorine  APT  as 
calculated  in  this  section. 

Further  Discussion  and  Interpretation  of  the  HOCi  APTs 
We  have  chosen  the  set  of  coordinates  shown  in  Fig.  4-1  to 
analyze  the  total  hydrogen  atari  APT  in  terms  of  the  CCFO  contrib- 
utions for  the  series  of  triatomic  molecules  studied  in  this  work. 
In  this  bond  coordinate  system  the  z axis  is  placed  along  the  OH 
bond  and  the  x axis  is  directed  toward  the  perturbing  x atom.  The 
charge  tensor  is  diagonal  and  obviously  invariant  to  the  coordinate 
system  used.  The  hydrogen  atom  APTs  in  the  original  coordinates 
were  given  in  Tables  3-8  , 3-22  and  3-40.  The  APTs  in  the  bond 
coordinate  system  are  presented  in  Tables  4-3  to  4-6  arranged  to 

show  each  separate  contribution  to  P for  all  four  HOX  molecules , 

-x 

for  purposes  of  comparison.  Table  4-3  shows  that  the  main  changes 

in  the  total  hydrogen  atom  APTs  for  HOX  from  molecule  to  molecule 

are  observed  in  the  values  for  the  3p  /3x  and  3p  /3z  elements , 

x z 

as  might  be  expected  (see  Fig.  4-1) . Table  4-4  shows  that  there 

is  no  substantial  change  from  one  molecule  to  another  in  the  charge 

tensors  for  HOX.  However,  Tables  4-5  and  4-6  show  that  there  are 

substantial  differences  in  the  charge-flux  and  overlap  tensors 

particularly  in  the  3p  /3x(CF),  3p  /3z(CF),  and  3p  /3z(0)  elements. 

X z z 
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Fig.  4-1. 


Bond  coordinate  system  for  the  hydrogen  atom  of  the 
HOX  series  of  molecules,  X = H,  Ci,  0,  and  F. 
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TABLE  4-3 

AB  INITIO  4-31G  APTs  ON  THE  HYDROGEN  ATOM  FOR  HOX 

(UNITS  ARE  e) 


P (Total) 

X 

—x 

/ 0.34 

0 

-0.02  \ 

• 

0.47 

* 

H 

-0.10 


(0.21 

° 

-0.03 


( 


0.29 

0 

-0.04 

0.25 

0 

-0.03 


0 

0 

0.49 

0 

0 

0.47 

0 

0 

0.50 

0 


0.09 


c£ 


F 


Note:  The  coordinate  system  is  shown  in  Fig.  4-1. 
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TAB^E  4-4 

THE  CHARGE  CONTRIBUTION,  PH  (C)  TO  THE  HYDROGEN  ATOM  APT 
FOR  HOX  (SNITS  ARE  e) 


PH(C) 


X 


0 

0 

\ 0 


.41 


/ 0.44 
0 

\ ° 

/ 0.45 
0 

\ 0 


0 

0.41 

0 

0 

0.44 

0 

0 

0.44 

0 

0 

0.45 

0 


H 


c£ 
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TABLE  4-5 

THE  CHARGE-FLUX  CONTRIBUTION  PH(CF)  TO  THE  HYDROGEN  ATOM  APT 

FOR  HOX  (UNIX'S  ARE  e) 


H 

P (CF)  X 

— X 


( -0 


.08 


\ -°* 


06 


0 

0 

0 


(-0.24 

° 

0.05 


0 

0 

0 


-0.06 

0 

0.12 


(-0.12 

° 

0.04 


0 

0 

0 


(-0.17 

° 

0.05 


0 

0 

0 


H 


Cl 


0 


Note:  The  coordinate  system  is  shown  in  Fig.  4-1. 
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TABLE  4-6 

THE  OVERLAP  CONTRIBUTION  PH (O)  TO  THE  HYDROGEN  ATOM  APT  FOR  HOX 

(uni$s  ARE  e) 


PH(0) 
— x 


X 


(0.01 
0 

-0.04 


(0.004 

° 

-0.08 


(-0.03 

° 

-0.08 


(-0.03 

° 

-0.08 


0 

0.06 

0 

0 

0.04 

0 

0 

0.03 

0 

0 

0.04 

0 


-0.03  \ 

0 

-0.38  j 


H 


Cil 


O 


F 


Note:  The  coordinate  system  is  shown  in  Fig.  4-1. 
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Let  us  first  consider  the  overlap  tensor  elements  of  the 

hydrogen  atom  of  the  HOX  series,  shown  in  Table  4-6.  A more 

careful  inspection  of  this  Table  shows  that  the  diagonal  element 

9p  /9z  appears  to  be  correlated  with  the  Mulliken  electronegativity 
z 

of  the  X atom.  This  correlation  appears  to  be  valid  even  for 

the  overlap  tensor  element  9p  /9z (0)  for  the  hydrogen  atom  of  HOC#.. 

z 

This  correlation  plot  is  shown  in  Fig.  4-2. 

Consider  the  reasonning  that  has  led  to  this  plot.  First  we 
observed  that,  with  the  exception  of  the  diagonal  9pz/9z(0)  elements 
of  the  hydrogen  atom  for  this  series  of  HOX  molecules,  all  the  other 
elements  of  the  overlap  tensors  are  quite  small  and  approximately 
constant  from  one  molecule  to  another.  However  the  diagonal  9pz/3z 
tensor  elements  change  from  -0.28e  for  the  hydrogen  atom  of  HOH  to 
-0.38e  for  the  hydrogen  atom  of  HOF.  We  observed  that  if  we  plot 
each  zz  element  of  the  overlap  tensor  versus  the  electronegativity 
of  the  x atom  we  obtain  a beautiful  straight  line  shown  in  Fig.  4-2. 

The  corresponding  overlap  tensor  element  for  the  hydrogen  atom  of 
CH^OH  is  also  included  in  that  figure  in  order  to  show  that  this 
element  from  this  molecule  also  fits  the  electronegativity  correlation. 
This  linear  electronegativity  correlation  was  also  observed  for  all 
of  the  elements  of  the  charge-flux  tensors  of  the  hydrogen  atom  of 
the  HOX  series  of  molecules  (shown  in  Table  4-5)  with  the  exception 
of  the  charge-flux  tensor  elements  for  the  hydrogen  atom  of  HOC#-. 

These  additional  electronegativity  correlations  are  shown  in  Fig.  4-3 
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We  therefore  used  the  electronegativity  correlation  of 
the  charge-flux  tensor  elements  of  the  hydrogen  atoms  shown  in 
Fig.  4-3  to  4-6  to  obtain  the  "predicted"  charge-flux  tensor 
for  the  hydrogen  atom  of  HOC#,  by  reading  the  charge-flux  tensor 
elements  directly  from  the  electronegativity  of  the  chlorine  atom. 
This  predicted  tensor  was  added  to  the  charge  and  overlap  tensors 
as  calculated  with  the  4-31G  basis  set  to  obtain  the  total  APT 
for  the  hydrogen  atom  of  HOC#,.  The  comparison  between  the  4-31G 
and  this  modified  APT  for  the  hydrogen  atom  of  HOC#-  is  shown  in 
Table  4-7. 

Next  we  rotated  the  coordinate  system  of  the  hydrogen  atom 
modified  tensor  to  the  coordinate  system  shown  in  Fig.  3-11. 

This  tensor  and  the  chlorine  APT  frcm  OC#,  given  in  Table  4-2  were 
used  to  obtain  the  oxygen  atom  APT  for  HOC#,  frcm  the  null  condition 
given  by  Eq.  (62).  The  resulting  APTs  for  HOC#-  are  shown  in  Table 
4-8.  Using  these  APTs  we  then  predicted  the  intensities  for  HOC#, 
as  listed  in  Table  4-9.  We  have  also  included  in  this  Table 
the  predicted  intensities  calculated  from  a 4-31G  basis  set  and 
from  a TZPP  basis  [76]  for  purposes  of  comparison.  This  Table 
shows  that  the  predicted  intensities  using  the  APTs  shown  in  Table 
4-8  are  in  remarkable  agreement  with  the  experimental  intensities. 

The  computer-simulated  spectrum  predicted  frcm  intensities  calculated 
from  the  modified  APTs  frcm  Table  4-8  is  compared  with  the  computer 
simulation  of  the  experimental  spectrum  of  HOC#,  in  Fig.  4-7. 
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TABLE  4-7 


H 


COMPARISON  BETWEEN  THE  4-31G  AND  MODIFIED  P FOR  HOC£ 


(UNITS  ARE  e) 


—x 


PH 

—x 


4-31G 


Modified 


P (C) 
—x 


P (CF) 
-x 


P (O) 
—x 


P (T) 
—x 


0 

0.44 

0 

0 

0 

0 


0.004  0 

0 0.04 

-0.08  0 


0 

0.49 

0 


/ 0.004 
0 

\ -0.08 


0 

0.44 

0 

0 

0 

0 

0 

0.04 

0 

0 

0.48 

0 


Note:  The  coordinate  system  is  shown  in  Fig.  4-1. 
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TABLE  4-8 

MODIFIED  APTs  FOR  HOC£  USED  TO  CALCULATE  PREDICTED  SPECTRUM 
IN  FIG.  4-7 (UNITS  ARE  e) 


P 


0 


— x 


(-0.64  0 

0 -0.32 

0 -0.07 


P 


H 


— x 


0 

0.19 

-0.06 


PC* 
— x 


0 

0.13 

0.13 


Note:  P°  Obtained  from  the  null  condition,  given  by  Eq.  (62) 
—x 

PH  (Modified) , see  text 
—x 

Pc  (Ab  initio  4-31G) , this  work. 

—x 

The  coordinate  system  and  atom  numbering  are  shown  in 
Fig.  3-11. 
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TABLE  4-9 

PREDICTED  INTENSITIES  FROM  AB  INITIO  AND 
MODIFIED  APTs  FOR  HOC5, 

(UNITS  ARE  km  mole  X) 


Band 

V (cm  ) 

A (4-31G) 

A (TZPP ) 

A (4-31G) * 

A (Expl . ) 

V (OH) 

3609.2 

75.3 

122.3 

41.9 

56.6 

V (A0) 

1239.9 

39.3 

45.8 

99.8 

73.9 

V(OCA) 

725.0 

0.089 

6.2 

9.5 

10.6 

Note:  See  also  Table  3-39.  A(4-31G)*  are  the  modified  4-31G 

intensities . 


Fig.  4-7.  Computer  simulation  of  the  predicted  spectrum  of 
HOC£  from  modified  4-31G  APTs  and  the  computer- 
simulated  gas  phase  spectrum  of  HOCJl. 

(a)  Predicted 


(b)  Experimental 


TRANSMITTANCE  TRANSMITTANCE 

■00  0.25  0.50  0.75  1.00  JJ.00  0.25  0.50  0.75  1.00 
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— I 1- 

357.50 


-+■ 


-+- 


H 1 1 1 1 1 1 1. 

315.00  272.50  230.00  187.50  U5.00 

FREQUENCY  (CM  (-1)  ) (X101  ) 


^00.00 


-1 1 1 

102.50  61 
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In  conclusion,  we  believe  we  have  demonstrated  how  power- 
ful these  techniques  for  predicting  spectra  can  be  in  the 
interpretation  of  the  spectra  of  free  radical  species  trapped 
in  inert  matrices  or  for  any  unstable  species  in  gas  phase. 

Valuable  predictions  can  be  made  either  by  using  APTs 
transferred  from  chemically  similar  molecules  or  by  using  APTs 
from  quantum  mechanical  calculations  made  by  ab  initio  SCF 
methods.  We  have  also  pointed  out  some  qualitative  and 
quantitative  applications  of  this  theory.  Most  of  our  absolute 
intensities,  either  quantum  mechanically  calculated  or  predicted 
by  proper  transfer  of  APTs,  are  within  factor  of  two  agreement 
with  the  experimental  values  and  therefore  often  within  the 
experimental  uncertainties.  For  free  radical  species  trapped 
in  a matrix  or  for  unstable  intermediate  species, it  is  more 
than  optimistic  to  expect  experimental  accuracy  to  be  better 
than  a factor  of  two  due  to  the  difficulties  in  handling  and 
measuring  the  integrated  absorbances. 

We  have  also  pointed  out  that  not  all  APTs  are  transferable. 

We  have  analyzed  the  changes  in  these  APTs  in  terms  of  the 
quantum  mechanical  CCFO  model  of  the  APT.  We  found  that  these 
changes  are  reflected  in  the  charge,  charge-flux  and  overlap 
contributions  to  the  APT  and  these  changes  often  reflect  on  changes 
in  the  chemical  environment  in  the  new  molecule  and  in  the  molecule 
to  which  the  APT  is  being  transferred.  We  have  also  pointed  out 
the  importance  of  an  additional  criterion  for  choosing  a quality 
basis  set  for  dipole  moment  derivative  calculations. 
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Finally,  the  success  of  our  calculations  indicates 
certainly  that  this  new  technique  can  be  one  of  the  most 
important  and  useful  techniques  for  the  qualitative  and 
quantitative  studies  of  free  radicals  or  unstable  intermediate 
species  trapped  in  inert  matrices.  It  has  important  analytical 
applications  in  the  study  of  the  photochemistry  of  those 
species,  as  well  as  in  studies  of  the  concentrations  of  these 
unstable  species  in  planetary  atmospheres. 
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